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Abstract 


Multivariable systems exhibit a higher degree of complexity than a single-variable 
system due to the strong coupling between the input and output variables, which makes the 
application of multivariable system control theory techniques difficult or impossible. On the 
other hand, the concept of fractional calculus has been successfully applied for the control of 
monovariable systems. This encourages researchers to consider introducing this concept to 
large systems of multiple inputs and outputs. However, research in this area is still at a very 
early stage. In this context, the main proposals of the thesis revolve around the following 
elements: 

- We propose a technique for the synthesis of fractional order controllers based on the 
optimization algorithms, in particular the Particle Swarm Optimization (PSO) algorithm for 
square multivariate systems with simplified decoupling with the aim of eliminate interactions 
so that inputs and outputs are controlled independently. 

-We propose two methods of analytical synthesis of fractional order controllers: the 
first method is based on the internal model control (IMC). The second method based on the 
Smith predictor (SP) structure. 

-In order to impose a fractional behavior on these systems, we use Bode's ideal 
function as a closed-loop reference transfer model for each independent loop. The command 
structure is equivalent to that of SP. It is also equal to the causality filter transfer function in the 
case of IMC-based control. 

For high order decoupled systems, they are approximated to type fractional order 
models (NIOPTD-I and NIOPTD-IT). We use the technique of direct synthesis to determine the 
parameters of the fractional controllers. The results obtained show the ability of the proposed 
control techniques to strengthen the control system, to eliminate the delay effect. They provide 


better robustness to noise and disturbances compared to other conventional methods. 


Keywords: Fractional order operator, fractional control system, multivariable systems, 


decoupling, modeling, synthesis. 


Résumé : 


Les systèmes multivariables présente un degré de complexité plus élevé qu'un système 
monovariables en raison du fort couplage entre les variables d'entrée et de sortie, ce qui rend 
l'application des techniques de la théorie de commande des systèmes multivariables difficile ou 
impossible. D'autre part, le concept de calcul fractionnaire a été appliqué avec succès pour la 
commande des systèmes monovariable. Ce qui encourage les chercheurs à penser d'introduire 
ce concept aux larges systèmes de plusieurs entrées et de sortie. Cependant, la recherche dans 
ce domaine est encore à un stade très précoce. Dans ce contexte, les principales propositions de 
la thèse s'articulent autour des éléments suivants : 

- Nous proposons une technique de synthèse des contrôleurs d'ordre fractionnaire basée 
sur les algorithmes d'optimisation, en particulier l'algorithme d'optimisation d'essaim de 
particules (PSO) pour les systèmes multivariables carrés avec découplage simplifiée dans le but 
d'éliminer les interactions afin que les entrées et les sorties soient contrôlées indépendamment. 

-Nous proposons deux méthodes de synthèse analytique des contrôleurs d'ordre 
fractionnaire: la première méthode est basée sur la commande à modèle interne (IMC). La 
deuxième méthode basée sur la structure de prédicteur de Smith (SP). 

-Afin d'imposer un comportement fractionnaire à ces systèmes, nous utilisons la 
fonction idéale de Bode comme modèle de transfert de référence en boucle fermée pour chaque 
boucle indépendante. La structure de commande est équivalente à celle de SP. Elle est 
également égale à la fonction de transfert de filtre de causalité dans le cas d'une commande 
basée sur IMC. 

Pour les systèmes découplés d'ordre élevé, ils sont approximé aux modèles d'ordre 
fractionnaires de type (NIOPTD-I et NIOPTD-I). Nous utilisons la technique de synthèse 
directe pour déterminer les paramètres des contrôleurs fractionnaires. Les résultats obtenus 
montrent la capacité des techniques de commande proposées à renforcer le système de 
commande, à éliminer l’effet de retard. Ils assurent une meilleure robustesse au bruit et aux 


perturbations en comparaison par les autres méthodes classique. 


Mots-clés : Operateur d'ordre fractionnaire, système de commande fractionnaire, systèmes 


multivariable, découplage, modélisation, synthèse. 
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Introduction générale: 


Les systèmes  multivariable (système à plusieurs entrées et plusieurs sorties 
interconnectées entre elles) présente un degré de complexité plus élevé qu'un système 
monovariable (système à une seule entrée et une seule sortie) ce qui rend l'application des 
techniques de la théorie de commande des systèmes multivariable difficile ou imopossible à 
cause de l'interconnexion entre les entrées et les sorties. 

Le principal défi pour les chercheurs est d'améliorer les performances du contrôle des 
processus MIMO et, si possible, d'éliminer l'impact de ces interactions. Parmi les solutions 
proposées pour résoudre ces problèmes, il est nécessaire l'analyse des interactions comme 
première étape [1]. Cette démarche permet de quantifier ou bien évaluer le degré des 
interactions dans un système multivariable, affin d'exprimer les difficultés créées par les 
éléments du couplage. Différentes méthodes d'analyse des interactions ont été développées 
[2] qui nous donne la meilleure configuration de commande où les interactions sont moins 
importantes. 

Dans la littérature, la méthode des gains relatifs (Relative Gain Array, RGA) proposée 
par Bristol en 1966 [1] est considérée comme une référence de base pour évaluer le niveau 
d'interaction entre les boucles des systèmes MIMO, Vu la simplicité de calcul, cette méthode 
est souvent utilisée dans l'industrie en association avec l'indice de Niederlinski [3 |. 

A l'étape d'analyse de l’interaction, si le niveau d'interaction est négligeable, le 
problème est résolu. Mais si le niveau d'interaction est élevé ou indésirable, dans ce cas les 
interactions dans le système doivent être affaiblies en réalisant un ensemble de synthèses 
mono variables en introduisant un découpeur [4-5-6] entre le processus qui doit être contrôlé 
et le contrôleur. De sorte que les entrées et les sorties soient contrôlées indépendamment. 
Cette dernière technique de commande est la plus utilisée en industrie, vu la simplicité 
d’implémentation et de l’entretien, ainsi que les performances qu’on peut atteindre en 
s’intéressant aux performances des boucles individuelles. Cette technique est appelée la 
commande décentralisée ou la commande multiboucles [6-7]. 

D'un autre côté et du point de vue commande, les contrôleurs PID classiques ont fait 
des progrès acceptables dans le contrôle des processus industriels. Ils sont considérés comme 
l'option la plus populaire et la plus utilisée en raison de leur simplicité et de leur facilité 
d'utilisation. Cependant, les contrôleurs PID ne peuvent pas toujours atteindre les 
spécifications souhaitées, Ces derniers temps, l'utilisation du calcul fractionnaire s'est 


développée dans plusieurs domaines par les chercheurs a fait apparaitre deux nouveaux 


concepts dans la théorie de la commande des systèmes à savoir : la commande d'ordre 
fractionnaire (COF) et les systèmes d'ordre fractionnaire (SOF). 

En 1990, Podlubny [8] a introduit de nouveaux contrôleurs PID d'ordre 
fractionnaire noté PIY"D# qui incluent des opérateurs de dérivation et d’intégration d'ordre 
fractionnaire sont une extension des contrôleurs PID classiques. 

Cependant, les méthodes de synthèse de lois de commande d'ordre fractionnaire 
s'avère très complexe car ils possèdent cinq paramètres à régler ou à calculer. Nombreux 
travaux de recherche et solutions ont été proposés pour développer des techniques de synthèse 
des contrôleurs fractionnaire, soit des méthodes optimales [9] basées sur les algorithmes 
génétiques et soit des méthodes analytiques pour la commande des systèmes monovariables 
d'ordre entier [10]. 

Récemment, Cette dernière méthodes de synthèse est la plus développé en industrie, 
en particulier le nouveau type de contrôleurs introduit par Bettayeb et Mansouri [11, 12], en 
combinant entre la commande à modèle interne (IMC) et la commande CRONE (Commande 
Robuste d'Ordre Non Entier) basée sur la fonction idéale de Bode imposée à la boucle fermée 
comme un modèle de référence (filtre de robustesse) pour la synthèse du contrôleur d’ordre 
fractionnaire en utilisant le principe de la commande IMC. Le nouveau contrôleur résultant de 
l'équivalence entre la structure de commande par modèle interne IMC et la structure classique 
à retour unitaire peut être décomposé en un contrôleur PID d'ordre fractionnaire cascadé avec 
un filtre d'ordre fractionnaire. 

De plus, dans le cas les systèmes monovariables contient un temps de retard important, 
différentes stratégies de contrôle appropriées sont utilisées. La meilleure d'entre elles est la 
stratégie de synthèse analytique de prédicteur de Smith (SP) a été adopté par Marzieh Safaei 
et Saeed Tavakoli [13], qui vise à combiner entre commande du prédicteur de Smith (SP) et 
la commande CRONE basée sur le choix du modèle de référence. Le nouveau contrôleur 
résultant peut être décomposé en un contrôleur PID d'ordre entier cascadé avec un filtre 
d'ordre fractionnaire. 

Le concept de calcul fractionnaire a été appliqué avec succès pour la commande des 
systèmes monovariable. Ce qui encourage les chercheurs à penser d'introduire ce concept aux 
larges systèmes de plusieurs entrées et de sortie. Cependant, en se référant à la littérature, 
nous constatons que la commande d'ordre fractionnaire des systèmes multivariables n'a pas 
été suffisamment discutée, car on peut dire que la recherche dans ce domaine est encore à un 
stade très précoce. 


Le travail présenté dans cette thèse s'inscrit dans le cadre du développement de stratégies 


appropriées pour la conception de contrôleurs d'ordre fractionnaires pour les systèmes 
multivariables carrés (les systèmes à plusieurs entrées et plusieurs sorties carrés.). L'objectif 
consiste à proposer des méthodes efficaces de synthèse du contrôleur d'ordre fractionnaire 
avec une facile implémentation et meilleur description de comportement du système en 
boucle fermée. Dans ce contexte, nous proposons dans notre contribution deux nouvelles 
méthodes de synthèse des contrôleurs d'ordre fractionnaires à savoir, les méthodes de 


synthèse optimale et les méthodes de synthèse analytique pour les systèmes multivariables. 


Cette thèse est décomposée en cinq chapitres organisés comme suit : 


-Le chapitre 1 introduit les notions mathématiques de base du calcul fractionnaire et 
des opérateurs d’ordre fractionnaires, et aussi la transformée de Laplace des dérivées et 
intégrales d'ordre fractionnaire. Nous allons aussi présenter quelques méthodes la réalisation 
des opérateurs d’ordre fractionnaire par les d’approximation numérique par exemple 


l'approximation basé sur la représentation diffusive, appuyée par des exemples de simulations. 


-Dans le chapitre 2, nous commençons par résoudre le problème des interactions pour 
les systèmes multivariables, après avoir analysé et quantifié ces interactions par la méthode 
RGA, nous présentées trois différentes technique de découplage pour affaiblir ou éliminer 
l’effet des interactions et sont: le découplage inversé, le découplage simplifié et le découplage 
idéal. Aïnsi, le nouveaux schéma de commande pour la commande des systèmes MIMO 
appelé la commande décentralisé et cela garantit les performances et les robustesses désirées 


car les entrées et les sorties sont contrôlées indépendamment. 


-Le chapitre 3 portera sur l'application de la méthode de synthèse optimale proposée 
pour les contrôleurs PI d'ordre fractionnaire (FO-PT) décentralisé par un algorithme génétique 
est implémenté sur un modèle réel de la colonne de distillation binaire [14-15]. La fonction 
de transfert non standard des opérateurs d'ordre fractionnaire est réalisée au moyen d'une 
approche diffusive. Cette approche est réalisée pour la première fois sur les systèmes 
multivariables. Une discussion des résultats obtenus révélera les améliorations avec celles 
obtenues grâce à l'utilisation du régulateur PID classique. Pour prouver l'efficacité de la 
méthode de synthèse proposée, une seconde simulation est réalisée sur un système MIMO de 
trois entrées et trois sorties (3x3). 


-Le chapitre 4 traite de la synthèse analytique proposé de commande fractionnaire pour les 


systemes multivariables de deux entrée et de deux sortie (systèmes TITO). Nous avons 
proposé une manière d'étendre la structure de commande par modèle interne IMC pour les 
systèmes SISO pour inclure à la commande des systèmes multivariables avec amélioration les 
performances du système en boucle fermée. L'effet de toutes les boucles est considéré comme 
une perturbation sur la boucle de commande considérée, d'où l'idée d'utiliser découplage , 
pour séparer les interactions entre les variables du système TITO et l’obtention les boucles 
indépendantes ( systèmes découplées de type SISO). tandis que l'étape suivante est une étape 
de modélisation dans laquelle les boucles indépendantes sont approximée avec des nouveaux 
modèles d'ordre fractionnaire a été proposé de type I et II [16-17] et est appelé: modèle 
d'ordre non entier plus un temps de retard de type I ou II (NIOPTD-I et NIOPTD-IT). Après 
ces étapes, les méthodes de synthèse du contrôleur d'ordre fractionnaire précédent basé sur la 
commande de IMC deviennent réalisables et appliquée à chaque boucle séparément. Nous 
utilisons également la la fonction idéale de Bode en boucle fermée dans la domaine 
fréquentiel comme modèle de référence dans la synthèse du contrôleur IMC-FOF-FOPI multi- 


boucle d'ordre fractionnaire. 


-Le chapitre 5, Englobera nos contributions le travail sur le développement de la 
synthèse des contrôleurs d'ordre fractionnaires à partir de modèles d'ordre entier avec retard 
vers des modèles d'ordre non entier à retard dans un schéma de Prédicteur de Smith. Dans ce 
contexte, nous proposons un modèle d'ordre fractionnaire de type II qui est équivalent aux les 
systèmes entiers d'ordre élevé comme première étape, dans notre prochaine étape nous 
utilisons la structure de la fonction de transfert idéal de Bode en boucle fermée comme 
structure équivalente au schéma de prédicteur de smith. Le nouveau contrôleur résultant peut 
être décomposé en un contrôleur PID d'ordre fractionnaire cascadé avec un filtre d'ordre 
fractionnaire. Nous appliquons la approche de la synthèse directe qui permet d'obtenir les 
paramètres du contrôleur FOPID directement en utilisant les paramètres des modèles d'ordre 
fractionnaire. et les paramètres du filtre d'ordre fractionnaire sont facilement obtenus en 
utilisant les formules des spécifications désiré du domaine temporel de la fonction de transfert 
idéal de Bode en utilisé comme modèle de référence. A la fin du chapitre, le schéma de 
commande proposé est implémenté dans les systèmes monovaribales avec un retard et sera 
généralisé dans les systèmes complexes (MIMO) de deux entrée et de deux sortie à retard, 
pour lequel nous montrerons l'efficacité du notre contribution, et la robustesse des lois de 
commande proposées qui élimineront l'effet du retard aux les systèmes commandé. 


L'approche de synthèse que nous proposons a été réalisée avec succès pour la première fois 


dans les systèmes monvariables ou multivariables. 


Finalement, nous terminerons notre travail par une conclusion générale et des perspectives 


potentielles de travaux futurs seront présentées à la fin du manuscrit de thèse. 
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Chapitre I : Calcul fractionnaire et équation différentielle. 


1.1 Introduction : 

Les systèmes d'ordre fractionnaire sont généralement décrits par des équations 
différentielles d’ordre non entier, qui modélisent réellement le comportement des systèmes 
physiques car ils donnent des modèles plus précis. Malgré l'énorme quantité de résultats 
publiés sur les équations différentielles fractionnaires et les systèmes dynamiques, nous 
pensons qu'il reste encore de nombreux problèmes ouverts. En effet, la théorie et l'application 
de ces systèmes sont encore des domaines de recherche très actifs car ils ont une longue 
mémoire (comportement héréditaire). 

L'objectif principal de ce chapitre est de présenter certaines bases théoriques des calcul 
fractionnaire pour le développement de la commande des systèmes dynamiques , tout en 
rappelant les principales définitions et propriétés des opérateurs d'ordre fractionnaire et leurs 


approximations qui existent. 


1.2 Calcul fractionnaire : 


Le calcul fractionnaire (CF) est une variante du calcul différentiel dans laquelle les intégrales et les 
dérivées sont d'ordre fractionnaire. Il peut être mis en œuvre pour des systèmes dynamiques qui ont été 
modélisés à la fois par des équations d'ordre entier et d'ordre fractionnaire. Le calcul fractionnaire a 
suscité un intérêt considérable dans la recherche récente et a trouvé diverses applications, notamment 
dans les systèmes de contrôle. La modélisation d'un système en ordre fractionnaire permet d'obtenir de 
meilleures performances par rapport à son modèle en nombre entier. Le modèle d'ordre entier est lui- 
même une approximation d'un modèle FO équivalent. La mise en œuvre d'un modèle fractionnaire est 
comparativement complexe, mais elle permet d'améliorer les performances telles que la robustesse, la 
stabilité transitoire, le filtrage du bruit et le rejet des perturbations. La conception d'un contrôleur pour 
un système donné est un autre problème majeur, et les contrôleurs d'ordre fractionnaire (FO) ont 
démontré des performances louables même dans ce domaine. Les contrôleurs PID d'ordre fractionnaire 
(FOPID) ont été mentionnés pour la première fois dans les travaux de Podlubny [18, 19,20] dans 


lesquels un intégrateur et un différentiateur FO sont conçus pour un système d'ordre fractionnaire. 


1.2.1 L'opérateur intégro-différentiel d'ordre fractionnaire 
Le calcul fractionnaire est une généralisation de l'intégration et de la différentiation à 
l'opérateur fondamental d'ordre fractionnaire aDŸ où a et t sont des limites de l'opération. 


L'opérateur intégro-différentiel continu est défini par : [20] 


net > 0 
at =: = (1.1) 
f (dr) a < 0 


où a € R est l’ordre de l’opération. 

Il existe de nombreuses définitions du calcul fractionnaire. Les trois définitions suivantes 
sont largement utilisées dans le domaine des systèmes de contrôle [21]. 
1.2.1.1 Définition de Grunwald-Letnikov (G-L) 


La définition de Grunwald-Letnikov (G-L) est exprimée comme suit: 


DEF = lime à D (7) it) (1.2) 


où, En est la partie entière et a et t sont les limites de l'opérateur , n est le nombre entier 


qui satisfait à la condition n—1<a<n et la valeur du coefficient binomiaux est donnée 


par: 
a (a +1) 
( ) — 3 (15) 
j TG +1)7(a-—j+1) 
La fonction Gamma utilisée dans l'équation ci-dessus peut être définie par l'équation suivante: 
ré l te Ted RC) > D (1.4) 
0 


La définition GL (Grunwald-Letnikov) est utilisée pour les évaluations numériques. Cette 
définition est très bénéfique pour obtenir une solution numérique des équations différentielles 
fractionnaires [22]. 

1.2.1.2 Définition de Caputo 

La définition de Caputo est surtout populaire dans les applications d'ingénierie [23], car elle 
est directement liée au type des conditions initiales et au type de la dérivée fractionnaire. 
Comme l'indique I. Podlubny [18], cette définition autorise des conditions initiales telles que 
y(0), y'(0), et non des conditions fractionnaires telles que y°((0). De plus, la dérivation 
fractionnaire de Caputo de la constante est bornée. Elle est définie par 

if fo 
& = 


D? = 
Hi T(n-a)]J (t-T)"+1 
to 


(1.5) 


où, n est un entier qui satisfait la condition n —1 <a <n, « est un nombre réel et 
t, et t sont les limites d'intégration. 


La définition de RL sert une notation générale de la dérivée générale, tandis que la dérivée de 


Caputo est une notation générale de la différenciation de manière classique [19]. 

1.2.1.3 Définition de Riemann-Liouville (R-L) 

Liouville a défini la dérivée d'ordre arbitraire comme une série infinie. L'inconvénient est que 
l'ordre doit être limité aux seules valeurs pour lesquelles la série converge. Il a appliqué cette 
théorie à la théorie du potentiel, et a été le premier à tenter de résoudre des équations 
différentielles par des opérateurs fractionnaires. Riemann a utilisé la généralisation d'une série 
de Taylor pour dériver une formule relative à l'intégration d'un nombre arbitraire. 

Il est défini comme suit 


7) [. Fée 
en —( dt GŒ=T)a-nti 


où, n est un entier qui satisfait la condition n—1<a<n ,a est un nombre réel, 1 est 


4 DE = D'INa f(E) = dt (1.6) 


l'opérateur intégral et t, et t sont les limites d'intégration. 
La définition de RL pour l'intégrale et la dérivée fractionnaire est appropriée pour trouver la 
solution analytique de fonctions simples comme et, t?, cos(t) [24]. 
1.2.2 Transformations de Laplace de la dérivée d'ordre fractionnaire 
La transformée de Laplace de la fonction f(t) est la fonction F(s) de la variable complexe s 
définie par : 
F(s) = LE] = Je stf(bdt (17), 

Pour l'existence de l'intégrale (1.7), nous avons besoin de définir la fonction d'ordre 
exponentielle. 
1.2.2.1 Définition 1: 
Une fonction f est d'ordre exponentiel a s'il existe des constantes M > 0 de sorte que pour 
certains T > O, 

[f()I<Met,t>T (1.8) 
Pour l'existence de l'intégrale (1.7), on a ce théorème. 
1.2.2.2 Définition 2: 
Si f est continue par morceaux sur [c; 1) et d'ordre exponentiel a, alors la transformée de 
Laplace L(f(t)) existe pour Re(s) > a et converge absolument. 
Afin d'appliquer la transformée de Laplace à des problèmes physiques, il est nécessaire 
d'invoquer la transformée inverse. si L(f(t)) = F(s), alors la transformée de Laplace 
inverse est notée par 

F[F(S)]=/f0, t30 (1.9) 


ui renvoie la transformée de Laplace d'une fonction à la fonction d'origine et est définie par: 
P 8 P 


foie fG)l= I estF(s)ds, c = Re(s) > © (1.10) 


où Co se trouve dans le demi-plan droit de la convergence absolue de l'intégrale de Laplace. 


1.2.3 Propriétés de la transformée de Laplace 

Nous mentionnons ci-dessus certaines propriétés de la transformée de Laplace dont nous 
aurons besoin plus tard. 

- La transformée de Laplace est une cartographie linéaire, c'est-à-dire pour toutes les fonctions 


f et g admettant des transformées de Laplace et pour tous les réels « et B : 


L(af() + Bg()) = aL(f &)) + BL(gG)) (1:11) 
- Sous les mêmes conditions, nous avons 
L(CF * g)@)) = L(F ())L(g(E)) (1:12) 


où * désigne le produit de convolution défini par 


t 


Grp = fr@-0 tar = | at -0 rer (1.13) 
0 


0 


avec f et g sont des fonctions continues sur [0, co). 
- Une autre propriété utile dont nous avons besoin est la transformée de Laplace de la dérivée 


d'un nombre entier d'ordre n d'une fonction f (€) : 


n-1 n=1 
L(F(E)) —_ s"F(s) _ st k-1F 0) (0) _ s'F(s) . sk fR-k-n (0) (1.14) 


où F(s) = L(f(E)). 
Pour l'instant, nous allons donner la transformée de Laplace de la dérivée fractionnaire d'ordre 
a > 0 vue précédemment. 


- L'écriture de l'intégrale fractionnaire d'ordre «4 > 0 (1.1) comme produit de convolution: 
pa-1 
IF) = Fo 10 (1.15) 
- Définition de Riemann-Liouville 
L'expression de la transformée de Laplace de la dérivée fractionnelle de Riemann-Liouville de 
d'ordre a > 0 est donnée par: 
L(D2()) = SF) 2 sep EG]. (1.16) 


Ou s représente la variable de Laplaceetn—-1<a<n 


- Définition de Caputo 


L'DPFOI=S FCO ESEESEEt7E 0), n—1<a<n (1.17) 
- Définition de Grunwald-Letnikov (G-L) 
L[ .DEFE)] = SF(S) (1.18) 


En utilisant la formule de la transformée de Laplace de la convolution (1.12) on obtient la 
transformée de Laplace de l’intégrale de Riemann-Liouville et celle de Gründwald-Leitnikov : 

LA EF) = À, a > 0} (1.19) 
Ou s représente la variable de Laplaceetn—-1<a<n 


1.3 Applications du calcul fractionnaire 


La CF est applicable dans diverses applications, dans presque tous les domaines de la science 


et de l'ingénierie, comme l'illustre la figure 1.1. 


Application 


1 | Médicale 
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=] 
= 
=] 
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© 
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A Conception et analyse de circuits | Ÿ 


5 | Analyse des systèmes dynamiques 


Figure 1-1:Applications du calcul fractionnaire 


Par exemple, les applications de CF sur les systèmes dynamiques peuvent être classées 


comme illustré dans la Figure 1.2 


Modlation 


Modélisation d'un système dvn se du modèle d'ordre Conception différents control implémentation pratique del: 
3 A! le 
modèle é fractionnaire urs pour 1 controleur concue 


[ractionnaire 


Figure 1-2:Applications de FC sur les Systèmes Dynamiques 


1.4 Approximation des opérateurs d’ordre fractionnaires 
La simulation d’un Opérateurs d’ordre fractionnaire " s’ et s "est dans la plupart des cas très 
compliquée [25-26]. Plusieurs techniques, permettant de simuler les opérateurs d’ordre fractionnaire, 
ont été développées dans la littérature. La plupart d'entre elles sont basées sur l'approximation dans le 
domaine "s” de la fonction irrationnel (d’ordre non entier) de l'opérateur par une fonction rational 
(d'ordre entier). Ces techniques consistent à calculer la sortie du système en utilisant un modèle 
rationnel continu équivalent à l’aide d’une représentation spéciale, [27]. Elles sont appelées des 
approximations analogiques ou des approximations du domaine fréquentiel [28]. Parmi ces méthodes 
on peut trouver : 
e La méthode de Carlson 
e La méthode de Matsuda 
e La méthode de Charef 
e La méthode EFC (Expansion Fractionnaire Continue) 
e La méthode d'Oustaloup 
e La Méthode diffusive 
Par conséquent, dans le troisième chapitre de cette thèse, nous choisirons la méthode diffusive 
pour approximer des opérateurs d'ordre fractionnaire pour la commande décentralisée proposée des 
systèmes multi variable, Mais dans les quatrième et cinquième chapitres, nous choisirons la méthode 
d'Oustaloup qui est plus largement utilisée dans l'industrie. 
1.4.1 Méthode d'approximation d’Oustaloup 
Cette méthode est basée sur l'approximation, en temps continu, de la fonction transfert irrationnel de 
l'opérateur d'ordre fractionnaire sŸ à une fonction transfert rationnel. Cette approximation utilise une 


distribution récursive de N zéros, et N poêles se trouvant dans la bande de fréquence [w, ,w}] , 


L’approximation par la méthode d’Oustaloup de l’opérateur fractionnaire sŸ est donnée par [29-30]: 


s+w; 


G(sS)=sh = sh, w,] =Kk] [= , HE R (1.20) 
L 


i=1 
avec : N = approximation d'ordre 
Où: le gain, les pôles et les Zéros peuvent être calculés respectivement à partir de : 


K = Gain = wi (62h) 


N+ (1H) 
ne 
w;! = Zéros = w, (2) és (122) 
Wp 
EN+ (14H) 
M SP à (7) FE (1.23) 
Wb 


1 
1.4.1.1 Approximation de l'intégrateur d'ordre fractionnaire; 


Nous avons choisis un ordre d’intégration À = 0.5 , et la bande de fréquence [w 


107%,w, = 10°], avec des changements dans le nombre de pôles (N = 4,N = 6etN 


20 ) Nous obtenons les résultats suivants : 


Amplitude(dB) 


Phase (deg) 


: 10 10< 10* 10 
Fréquences (rad/sec) 


10° 10* 10°* 10< 


Figure 1-3:Diagramme de Bode d'un intégrateur d'ordre fractionnaire approché par méthode 
d’Oustaloup. 


1.4.1.2 Approximation de dérivateur d'ordre fractionnaire s“ 


Nous avons choisis un ordre de dérivateur u = 0.5 et la bande de fréquence [wy = 


107 ,wn = 10°] , avec des changements dans le nombre de pôles (N = 4,N = 6etN = 


20 ) Nous obtenons les résultats suivants : 


Diagramme de Bode 


Ambplitude(dB) 


Phase (deg) 


Fréquences (rad/sec) | 


Figure 1-4: Diagramme de Bode d'un dérivateur d'ordre fractionnaire approché par méthode 
d’Oustaloup. 
On remarque que pour un bon choix de paramètres d'approximation on obtient une meilleure 
approximation d'opérateurs d'ordre fractionnaire, le gain et la phase sont très proches du gain et de la 
phase idéaux. 


1.4.2 Approximation avec la représentation diffusive 


Le principe fondamental de la représentation diffusive est permet d'exprimer une grande variété 
d'opérateurs intégro-différentielles d'ordre fractionnaire causaux de façon non héréditaire par des 
systèmes linéaires dynamiques entrée-sortie de nature diffusive. Etant donnée H(s) la fonction de 


transfert (non rationnelle) associée à l'opérateur causal de convolution H (d/ ae)a réalisation diffusive 


canonique de cet opérateur est exprimée, lorsqu'elle existe, par la réalisation d'état : 


D t) = —E(E, t) + u(t), V(E, t) — 0, Ë ER” 
y = JL uG, D V(E, D)dE 


Où u(®) est appelée représentation diffusive de H(s) 


(1.24) 


1.4.2.1 Représentation au moyen de la transformation de Laplace 
La représentation diffusive u(£&)d'un opérateur pseudo-différentiel à temps invariant de symbole 


H(s) est définie, lorsqu'elle existe [31-32-33-34-35-36], comme solution de l'équation intégrale : 
H(s) = f*" #O qe (1.25) 


0 s+8 


L'équation (1.25) est la fonction de transfert associée à la représentation diffusive. En d'autres 
termes, u(&) La représentation est obtenue directement par la transformation de Laplace inverse de la 
réponse impulsionnelle h(t) La représentation diffusive est donc une utilisation "à contre sens" de la 


transformée de Laplace dans laquelle t joue le rôle de la variable de Laplace : 


L£ £ 
Représentation diffusive u(£) + réponse impulsionnelle h(t) =. Symbole H(s) 


= Le : æ Et _1 
CE) + RCE) + H(s) ; Où: he) = L(UC)) = J, TS; u(E) = LT (RCE) (126) 
Par exemple le symbole diffusif de l'intégrateur fractionnaire de l'ordre y : 

â\ . 1 
H(2) où H()==z for 0<À<1 (1.27) 


Le symbole diffusif est exprimé par [37] : 
Ho = or £ > 0 où a est l'ordre d'intégration. 


Et la réalisation diffusive du contrôleur PIAD4 fractionnaire avec l'entrée uw et la sortie y peut 
s'exprimer par [9]: 
dpt) = —b(,t) +u(E) 
+00 


1.28 
70 = | rOUE,Dds Cas 


0 


La fonction de transfert du contrôleur PI}D4 fractionnaire au moyen d'une représentation diffusive 


est: 
ut 
C{s) = I d 1.29 
Mr (1:29) 
Où v la représentation diffusive sait comme suit: 
1 sin(An) sin(un) 
V(8) = KP CORRE : s) (1.30) 
i TT TT 
À 


1.4.2.2 Approximation numérique d’un intégrateur s7 
Nous avons choisis un ordre d’intégration À = 0.7 , la bande d’approximation [éin = 


1077, Emax = 107] , avec des changements dans le nombre de pôles (Q = 10 ,Q = 15 ,Q 


30 ), Nous obtenons les résultats suivants : 


Diagramme de Bode 


— — Réponse exacte 
: In L4 


Fe Lu = 


Amplitude(dB) 


Phase (deg) 


10710 10° 10° 10° 1019 
Fréquences (rad/sec) 


Entre 10 76-10 407660 =100=416;:0= 30) 
Figure 1-5:Diagramme de Bode d'un intégrateur d'ordre fractionnaire approché par 
représentation diffusive. 
1.4.2.3 Approximation numérique d’un dérivateur s“ 


Nous avons choisis un ordre dérivation uw = 0.7, et la bande d’approximation [&,,;, = 


1077, Emax = 107] , avec des changements dans le nombre de pôles (Q = 10 ,Q =15 ,Q 


30 ), Nous obtenons les résultats suivants : 
Diagramme de Bode 
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Figure 1-6: Diagramme de Bode d'un dérivateur d'ordre fractionnaire approché par 


représentation diffusive. 
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Les figures obtenues montrent que la représentation diffusive proposée assure une bonne 
approximation de l'intégrateur d'ordre fractionnaire car un bon choix de paramètres d'approximation 
garantie une bonne approximation dans un bande qui dépend des pulsations de système. On remarque 
que pour Emin = 1077, émax = 10”,u = 0.7 et Q = 30 on obtient une meilleure approximation de 
dérivateur fractionnaire d’ordre y = 0.7 le gain et la phase sont très proches du gain et de la phase 
idéaux (+14 dB et +63°). 
1.5 Modèles et représentations des systèmes d'ordre fractionnaire 
Le système d’ordre fractionnaire est généralement représenté par des modèles d’ordre non 
entiers. Il existe trois modèles d'ordre fractionnaire qui peuvent décrire le système dynamique 
en temps continu comme suit: Soit par l’équation différentielle d'ordre fractionnaire 
généralisée, ou, par représentation d’état fractionnaire ou bien, par fonction de transfert 
fractionnaire: 
1.5.1 Equation différentielle d'ordre fractionnaire 
Dans ce qui suit, nous nous concentrerons uniquement sur la représentation en temps continu. 
Les équations différentielles d'ordre fractionnaire d'un système dynamique en temps continu 
peuvent s'écrire comme suit 
HD.) (y, y) = G(DPOP:P2..Eu (uv, ……ux) (31) 
Avec : où y;,u; fonctions du temps représentatives des sorties et des entrées respectivement 
du système décrit par (1.31). 
H(), G()sont des combinaisons (pas nécessairement linéaires) de lois de l'opérateur dérivé 
d'ordre fractionnaire et a;,B; > O sont des dérivées d'ordre fractionnaire relatives 
respectivement à la sortie et à l'entrée du système. 
Pour les systèmes monovariables linéaires invariants SISO LTI FOS dans le temps, l’équation 
(1.31) devient : 

H(Dao%a2.aN)y(t) = G(DPoP1P2 su Bu }u(t) (1.32) 
On pose : 
{corne nn et ax ,bx ER (1.33) 
Finalement on obtient l’équation différentielle d’ordre fractionnaire : 
anDy(t) + an-1D M1 y(E) + + + aÿDy(t) = bn DPmu(t) + amn_1DPm-1u(E) + + 
boDou(t) (1.34) 
Selon la relation entre les ordres de dérivation fractionnaire de l’équation différentielle (1.34), 
on distingue deux types du système d’ordre fractionnaire: 


e Les systèmes d’ordre commensurables 


e Les systèmes d’ordre rationnels 
1.5.2 Fonction de transfert fractionnaire 
En appliquant la transformée de Laplace de l’équation (1.34) avec des conditions initiales 
nulles, les représentations entrée-sortie des systèmes d'ordre fractionnaire peuvent être 
obtenues. Dans le cas des modèles continus, un système d'ordre fractionnaire sera donné par 
une fonction de transfert de la forme: 


Y(s) _ DS PM + bym1SPM-1+..+p0s Po 
U(S) — Ansn+ansn+...+a9s40 


G(S) = (1.35) 


Comme on peut le voir dans les équations précédentes, un système d'ordre fractionnaire a une 
fonction de transfert d'ordre irrationnel dans le domaine de Laplace. Pour cette raison, on peut 
dire qu'un système d'ordre fractionnaire a une mémoire illimitée, et évidemment les systèmes 
d'ordre entier ne sont que des cas particuliers. 

Dans le cas d'un système d'ordre commensurable, la fonction de transfert en temps continu est 


donnée par 


n a\k 
G(s) = Dk=o be(s")" (1.36) 


20 ax (SE 
qui peut être considérée comme une fonction pseudo-rationnelle, H(1), de la variable À = 5°, 


D'où on peut considérer la fonction rationnelle H() du variable À = s® est : 


n k 
HA) = Zoe (137) 


7 JkeoakA 
Dans cette thèse, le comportement des systèmes fractionnaires étudiés sera abordé par des 
modèles commensurables avec les propriétés de linéarité et d'invariance en temps continu où 
la dérivée fractionnaire de Caputo est appropriée. 
1.5.3 Représentation d’état fractionnaire 
Pour les systèmes d'ordre fractionnaire commensurables, des critères puissants ont été 
proposés. En utilisant l'hypothèse l'ordre commensurable, le système de l’équation (1.36) 
admet également une représentation de type espace d'état. 
Dans ce contexte, Le modèle d’état d’un système d'ordre fractionnaire linéaire, continu 


invariant multivariable, s’écrit sous la forme : 


or = Ax(t) + Bu(t) (1.38) 
y(t) = Cx(t) + Du(t) 
Avec :avecu € R,y € R,a € R*; le vecteur d’ordre fractionnaire est : 
(a) = [& ea 4) (1.39) 
et 
Dx= [Dex Dax, … Dax] (1.40) 


Dans le cas particulier du système d’ordre commensurable à celle du système d’ordre entier 
par l’utilisation de l’opérateur D® et la définition : 
DOX = ii = 12.1 (1.41) 
En utilisant La transformée de Laplace de l'équation différentielle de l'équation (1.38), la 
fonction de transfert correspondante peut être exprimée sous la forme : 
G(s) = C[(sT1, — A)71]B (1.42) 

1.6 Fonction de transfert idéal de Bode 

Bode [38] a été le premier qui a proposé une fonction de transfert en boucle ouverte d'un 
système asservi à retour unitaire garantissant la propriété de l'iso-amortissement. Cette 


fonction connue sous le nom de fonction de Bode idéale est donnée par : 


Le (s) = (=) = BeR. (1.43) 


Te SB 
Où: tr. constante de temps, w. est la fréquence de coupure désirée et le paramètre £ la pente 
de la caractéristique idéale du gain. En fait, la fonction de transfert L(s) est un intégrateur 
d'ordre fractionnaire pour £ > 0 Les diagrammes de Bode de L, (s) (1 < B < 2) sont très 
simple (figure 1.8) La courbe d'amplitude est une ligne droite de pente 
constante - 20$ dB /dec, et la courbe de phase est un trait horizontal à —fBr/2 
rads. La courbe de Nyquist se réduit à une ligne droite passant par l'origine 
avec arg Lg (jw) = —-Bnr/2rad. 

Maintenant considérons le système a retour unitaire représenté sur la figure 1.7. 
avec Ls (s) est la fonction de transfert idéal de Bode. Ce choix de L; (s) donne un 
système en boucle fermé avec la propriété souhaitable d'être peu sensible à la variation du 


gain. 


E(S) 7 re (=) Y(S) 


Figure 1-7:Boucle de commande idéale de Bode. 


La fonction de transfert en boucle fermée : 


… Y(s) _ Lg (s) mn 1 …. 1 … B 
H(s) = no TES avec K = w. (1.44) 


(Pour : 1 < f < 2 en à un système d’ordre fractionnaire oscillatoire). 
Ce système à la propriété de robustesse, en termes d’invariance sous changement 


d’échelle de fréquence. 


20Log|L (o)| 


-206dB /dec 


OdB 


Log w 


PM =n(1-f /2) . 


Figure 1-8:Réponse fréquentielle de la fonction idéale de Bode en boucle 
Les caractéristiques générales de la fonction de transfert de Bode sont les suivantes [39]: 


a) En boucle ouverte : 
- Gain : une pente constante de — 20$dB/dec. 
- Fréquence de coupure : une fonction de K. 


- Phase : ligne horizontale de -B : 
- Nyquist : ligne droite avec un argument de - f = à 


b) En boucle fermée avec retour unitaire : 


- La marge du gain est infinie. 
- La marge de phase est constante :@, = B (1 = £) ; elle dépend seulement de f . 


-La valeur maximale de la fonction de sensibilité S(s). 
S(s) = max | = = 1<B<2 (1.45) 
1+L3 Gw) (©) 

Le comportement dynamique d'un système en boucle fermée dépend de son comportement 
en boucle ouverte. De ce fait, la boucle fermée résultante, dans le cadre d'une configuration 


typique à retour unitaire préserve la caractéristique de robustesse de la fonction idéale de 


Bode. Pour cette raison, ce système est considéré comme un modèle de référence par 
plusieurs stratégies de commande. Afin de pouvoir fixer les paramètres de ce modèle de 
référence en fonction de la dynamique désirée, une étude de la fonction idéale de Bode dans 
différents domaines est nécessaire. 

1.6.1 Modèle de référence dans le domaine temporel 


La réponse indicielle de la fonction idéale de Bode en boucle fermée est donnée par [40] : 


ei = — Z- )- i= Corn") = 1 Eg(-(west)") (1.46) 
Où L'T est la transformée de Laplace inverse et E4 (x) est la fonction Mittag-Leffler pour le 
cas 1 < PB < 2. 

La figure 1.5. montre les réponses indicielles de l'équation (1.46) pour w,, = 1 et différentes 
valeurs de B [(voir à [40-131]. Comme nous l'avons vu, la réponse indicielle devient plus 
oscillatoire lorsque le dépassement augmente de £. 


Réponse indicielle 
T T 


2 


B=1 


Temps (sec) 
Figure 1-9:Réponses indicielles de la fonction idéale de Bode en boucle fermée pour w,, = 1 


et différentes valeurs de f. 


D'autres auteurs tels Manabe [41] et Barbosa [40], ont proposé d'autres méthodes qui 
aboutissent à une expression donnant la valeur du dépassement, plus simple. Il suffit de tracer 
les valeurs numériques du dépassement pour différentes valeurs de B, comme le montre la 
figure (1.10). Ensuite la courbe ainsi obtenue est approximée par un polynôme d'ordre 2 en 
utilisant la fonction "cftool" de Matlab. Le polynôme obtenu est donné par l'équation 


(1.47 Javec une somme de l'erreur quadratique moyenne d'approximation (MSE) de 2 x 107. 


M,(%) = 73.9(82 — 1.6793B + 0.6756) 1<B<15 (1.47) 
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Figure 1-10: Dépassement Mp en fonction de f 
La même procédure est appliquée pour avoir l'expression du temps de réponse, en traçant 
ts (2%) X a et ts (5%) X w,, pour différentes valeurs de B(Figure 1.11), les deux courbes 


obtenues sont approximées et donnent lieu aux expressions suivantes [13]: 


0.7885P — 0.2693 ee 
B — 0.8673 PS 


t, (2%) *X we = (1.48) 


k ne 1.078 <B < 1.486 
B2-20128+1056 Fun 


0.8128 — 0.2036 


B — 0.8007 
71568 — 7.9 


B? — 1.303P + 0.2578 


ts (5%) X we = (1.49) 


ts Ug 


7 £#: (5%) 


Figure 1-11:Temps de réponse (t;) en fonction de f 


Remarque: 


> Si les performances désirées sont spécifiées dans le domaine temporel le paramètre £ 


5 


peut être déduit approximativement à partir du dépassement Mp désiré. 
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> La valeur approximative de Tc peut être déduite à partir des équations (1.48), 
(1.49) en fonction de la valeur de temps de réponse Tr désiré. 
> il a proposé que le contrôleur soit conçu de sorte que la fonction de transfert de la 


boucle de commande soit équivalente à la "fonction idéale de Bode». 


1.6.2 Modèle de référence dans le domaine fréquentiel 


Cependant, afin de faciliter la synthèse du contrôleur en utilisant les caractéristiques dans le 
domaine fréquentiel de la fonction idéale de Bode, nous utilisons une nouvelle écriture de la fonction 
idéale pour Bode est proposée par Batayeb, M et Mansouri [11-12], qui est la suivante: 


LAS : 1 


1 
Ls (s) = (= 0<y<1, APR OEM met 


= TeSY+1 ? 


La fonction de transfert en boucle fermée est : 
Y(s) _ Les) 1 1 


H 2 — "| — "|  —— —————  _— 
(s) E(s) 1+L8(s) 1+(5 41  1+7es*1 (1.51) 
We 
x 1 
Où We = Er 


Les caractéristiques fréquentielles de la fonction idéale de Bode sont obtenues en posant s = 
jw. 

Le module et l'argument sont respectivement donnés par : 

WC 

w 


Le Gw)lus : 20log = -20(y+ 1)log(—) (1.52) 


TQw }r*1 
; | 1 : TT 
P(Lr OW)),u : arg Er = —(y + D; (1.53) 


Le dépassement de sa réponse indicielle dépend du paramètre y , et la rapidité de la réponse 
dépend de la valeur de la constante de temps t.. Ces deux paramètres peuvent être fixés 
comme suit : 


La marge de phase ®,, est donnée par: 


1 


En remplaçant (Le Gw)),,, par son expression donnée par l'équation (1.52) et en tenant 


dm = + arg| 


compte de la relation (y + 1), on obtient : 
bn =T-Y+105 (155) 


La valeur de y peut être déduite facilement de l'équation (1.55), elle est donnée par : 


LT mn 
Len T/2 


La valeur de +. peut être déduite facilement de l'équation (1.50), elle est donnée par : 


= (1.56) 


1 


(1.57) 


Te 
Remarque: 
> Si les performances désirées sont spécifiées dans le domaine fréquentiel le paramètre 
y, peut être déduit approximativement à partir du marge de phase D, désiré. 
> La valeur approximative de T. peut être déduite à partir des équations (1.56), 
(1.57) en fonction de la fréquence de coupure wc désiré. 
> 


1.7 Les Contrôleurs d'ordre non entier (Contrôleur fractionnaire) 


Le régulateur PID, appelé aussi correcteur PID (proportionnel, intégral, dérivé) est utilisé 
depuis plusieurs décennies dans les industries pour des applications de contrôle de processus 
[42]. La raison de leur grande popularité réside dans la simplicité de la conception et les 
bonnes performances, notamment le faible pourcentage de dépassement et le petit temps de 
stabilisation. Le régulateur PID d'ordre fractionnaire est une extension du contrôleur PID 
classique. Il est moins sensible aux changements des paramètres d'un système commandé, ce 


qui donne une amélioration de la robustesse. 


1.7.1 Contrôleur PI d'ordre fractionnaire 


Il est de deux types, 


e PI 
K; 
CS) =K, + (1.58) 
L'autre forme de PI est 
K 
C(s)=K, (: + +) (1.59) 


(1.58) et (1.59) sont liés comme suit : K, = K1, K, = K;-K; 
Ces deux formes de contrôleur PI} sont utilisées dans cette thèse. 


° [PI} 


C(s) = (x, n ) (1.60) 


L'autre forme de [PI]? est : 
K 
C(S)=K, (: +2) 


d 1 
(1.60) et (1.61) sont liés comme suit : K, = K,2  , K, = K2:* K; 
Où, 12>0  ; KRSER, 


Pour À = 1, (1.58) et (1.60) représentent un contrôleur PI entier, 
1.7.2 Contrôleur PD d'ordre fractionnaire P D" 


Il est de deux types, 
e PD# 
C(s) = K, + Kas' 
L'autre forme de PDF est : 


C(s) = K,(1 + K,s#4) 


(1.62) et (1.63) sont liés comme suit : K,, = K, Ka = Ki * K2. 
° [PDF 
C(s)= (K, + Kas) il 
L'autre forme de PDF est : 
C(s) = K,(1 + K,s)" 
(1.64) et (1.65) sont liés comme suit: K,, = ke Ki = Ka * Ko. 


où,u > 0; K,, Ka, ER. 


Pour y = 1, (1.62) et (1.64) représentent un contrôleur PD entier, C(s) = K, + Kas 


1.7.3 Contrôleur PID d'ordre fractionnaire 


e PI D4 
Ki 
C(s)=K, + — + Kist 
(s) D sÀ as 

L'autre forme de PI}D4 est : 

K 

C(s) . K; (: + 4 + Kst) 

(1.66) et (1.67) liés comme suit: K,, = K1, Ki = Ki : K, et Ka = Ki * Ka. 


Ces deux formes de contrôleur PID4 sont utilisées dans cette thèse. 


Où, 1>0,u>0 : K Ki Kay ER. 
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(1.61) 


(1.62) 


(1.63) 


(1.64) 


(1.65) 


(1.66) 


(1.67) 


Pour À = 1 et y = 1, (1.66) représentent un contrôleur PID entier, C(s) = K, + “ + K;s 


Suivant la variation de ces deux paramètres (1 et 4), on peut distinguer différentes possibilités 


des contrôleurs d’ordre fractionnaires (figure.1.12). 


PD PID 


PI 


A=1 À 
Figure 1-12:Possibilités des contrôleurs à trois actions 
*Sik=letu=l, alors c'est un contrôleur PID classique. 
*SiÀA=0etu=]1, alors c'est un contrôleur PD classique. 
*SiA= let pu = 0, alors c'est un contrôleur PI classique. 
*SiA=0etp= 0, alors c'est un contrôleur P classique. 


La structure générale d'un contrôleur PID d’ordre fractionnaire est illustrée à la Figure. (1.13). 


E(p) 


Figure 1-13:Structure de contrôleur PI D4 fractionnaire. 


Dans le régulateur PID d'ordre fractionnaire (FOPID), les opérateurs d'intégration et de 
dérivation 1 et D sont généralement d'ordre fractionnaire. Par conséquent, en plus de réglage 
des coefficients K,,,K; et K; , nous devons également de réglage deux autres paramètres, à 
savoir l'ordre d'intégration fractionnaire À ainsi que l'ordre dérivée fractionnaire 1. 

La figure 1.14 illustre la structure générale d'un système en boucle fermée, dans lequel un 
contrôleur PID fractionnaire est utilisé. Les contrôleurs fractionnaires ont deux paramètres de 
plus que le contrôleur PID conventionnel, ainsi deux spécifications supplémentaires peuvent 


être satisfaites, ce qui améliore la performance du système global. On peut s'attendre à ce que 


le contrôleur PI} D# puisse améliorer la performance du système. 


Figure 1-14:Système de commande d'ordre fractionnaire en boucle fermée. 


1.7.4 Avantages des contrôleurs PID d'ordre fractionnaire 


Les caractéristiques d'un contrôleur PID fractionnel sont résumées ci-dessous : 


Cinq spécifications différentes peuvent être atteintes, ce qui n'est pas possible dans le 
cas du contrôleur PID classique ce qui donne une grande flexibilité. 

Le contrôleur PID fractionnaire peut facilement atteindre la propriété d'iso- 
amortissement, par rapport au contrôleur PID classique. 

Pour les systèmes d'ordre supérieur, les performances du contrôleur PID classique se 
détériorent, tandis que le contrôleur PID fractionnaire fournit de meilleurs résultats. 
Pour un système à long retard, un contrôleur PID fractionnaire fournit de meilleurs 
résultats que le contrôleur PID. 

Un contrôleur PID classique offre une stabilité moins robuste, alors que le contrôleur 
PID fractionnel est plus robuste et stable. 

Généralement, un système non linéaire est linéaire à différents points de 
fonctionnement, et les contrôleurs sont conçus pour différents points de 
fonctionnement, alors qu'un seul contrôleur PID fractionnaire est suffisant pour un 


système non linéaire [43]. 


1.8 Conclusion : 


Nous avons présenté dans ce chapitre quelques définitions de base pour l’intégration et la 


dérivation d’ordre fractionnaire, ainsi que sa représentation dans le domaine de Laplace. Nous avons 
: 


introduit le concept du système d'ordre fractionnaire et les contrôleurs PID d'ordre fractionnaire. Nous 


avons présenté différents méthodes d’approximation des opérateurs d’ordre fractionnaire (méthode 


d'Oustaloup et la représentation diffusive). La qualité de l’approximation d’opérateurs fractionnaires 


dépend du choix des paramètres de l’approximation numérique. 


Chapitre 2: Etude des systèmes multivariables. 


2.1 Introduction 


Le phénomène d'interactions (ou: couplages ) entre les variables d'entrée et de peut être 
considéré comme la principale caractéristique qui distingue les systèmes multivariables des 
systèmes monovariables sortie [44-45-46-47]. Ces interactions peuvent être vus comme l'effet 
de plusieurs entrées sur une même sortie, ou l'effet d'une entrée sur de plusieurs sorties. En 
effet, l’analyse des interactions est une première étape très importante pour évaluer le niveau 
des interactions dans un système multivariable [1-48]. Ensuite, différents structures de 
commande sont proposés soit par la généralisation des techniques des systèmes SISO [49-50] 
ou bien élaborer des schémas de commande spécifiques aux systèmes MIMO telles que la 
commande par découplage (ou: commande décentralisée) [51-52] et la commande 
multivariable (ou: commande centralisée). Le schéma de commande pour commander les 
systèmes MIMO vise à éliminer, si possible, où réduire l'influence de ces interactions sur les 
performances du système en boucle fermée, pour rendre les entrées et les sorties indépendante 
ce qui garantit les performances et les robustesses désirées. 

Dans ce chapitre, nous commençons par définir et illustrer le phénomène des interactions 
dans un système multivariable. Ensuite, les méthodes de découplage des systèmes 
multivariables les plus rencontrées dans la littérature (à deux/trois entrées et à deux/trois 
sorties (2x2) et (3x3)). Nous présentons les différentes approches de conception des 
correcteurs conventionnelle (PI/PID) des systèmes découplés lorsqu'ils sont contrôlés par la 
structure de commande multi-boucle (decentralized). Le chapitre se termine par un aperçu sur 


la synthèse du régulateur PI/PID décentralisés. 


2.2 Types et méthodes d’analyse des interactions : 


Les systèmes multivariables peuvent être considérés comme des systèmes avec plusieurs entrées (u1, 
U2,... Un) et plusieurs sorties (y1, V2, Ym) [53-54-47] comme le montre la figure (2.1). 

Les interactions ou couplages peuvent être vus comme l’influence d’une entrée (références, 
commandes ou perturbations) sur plusieurs sorties, ou bien l’influence de plusieurs entrées 


(références, commandes ou perturbations) sur une même sortie. 


Figure 2-1:Système multivariable 


Dans le cadre de cette thèse, seuls les systèmes carrés seront considérés, c'est-à-dire que le nombre 
d'entrées est égal au nombre de sorties (n = m = p). 
Pour un système multivariable à deux entrées et deux sorties (2x2), le schéma de commande 
multivariable est illustré par la figure 2.2. Le procédé G(s) est commandé par deux régulateurs 


monovariables C11 (sS)et C22(S). 


| 
| 
| 
| 
l 
L 


Figure 2-2:Commande Décentralisée 


L’analyse des interactions permet également d’estimer s’il est nécessaire de prendre en compte les 
interactions durant la synthèse des correcteurs monovariables (synthèse multiboucle) ou non (synthèse 
monoboucle). 

Les objectifs d’une analyse des interactions sont donc multiples : 

- Choix de la stratégie de commande. 

- Choix de la configuration commandes-sorties dans le cas de la commande décentralisée. 

- Choix de la stratégie de réglage des correcteurs dans le cas de la commande décentralisée. 

Pour un système multivariable qui comprend m boucles, on aura m ! configurations possibles dont une 
seule est prise en considération. Plusieurs méthodes d’analyse des interactions ont été développées 
permettant de choisir la meilleur connexion et/ ou de juger le niveau d’interaction entre les boucles de 


régulation [55-56-57]. 


2.2.1 Analyse des interactions par la méthode de RGA (Relative Gain Array) 


Bristol [1-58] a développé une technique organisée pour l'étude des problèmes d’interactions de 
boucles de contrôle du processus multivariables basée sur la théorie d'un gain relatif. Pour un 
processus MIMO avec « n » variables, le RGA ( d;; ) entre la variable contrôlée « y; » et la variable 
manipulée « u; » est défini comme le rapport de deux gains en régime permanent s'exprime comme 


suit : 


Yi ) 
à. ( Ou; \ _ le gainstatique boucle ouverte 12 21 
PET (OL ) ” le gain statique en boucle fermée NRA L en 
U; 
de: 


dans l'équation. (2.1), le symbole de You.) est la dérivée partielle qui est calculée avec toutes les 
J u 


variables manipulées sauf que u; est constant. Par conséquent, ce terme est le gain en boucle ouverte 


(ou gain statique) entre y; et u; , qui correspond aux éléments de la matrice de gain k;;. De même, 


( Yi u) est évalué avec toutes les variables de sortie (contrôlées) à l'exception de y; qui reste 
j 
y 


constant. Ceci peut être réalisé en pratique en ajustant les autres variables d'entrée (manipulées) en 


0yi 


utilisant d'une procédure de contrôle intégrée. En conséquence, ( u) peut être obtenu comme 
J/u 


un gain en boucle fermée indiquant l'effet de u; sur y; tandis que toutes les autres variables définies 
sont maintenues constantes. Il est pratique d'organiser le gain relatif sous forme de matrice également 


connu sous le nom de La matrice des gains relatifs (RGA) qui est donné par l'équation. (2.2). 


A1 2 ln 
A=| As 2 ln (2.2) 
An An2 . Ann 


2.2.2 Calcul du RGA pour un système (2x2) : 


pour un système (2x2), supposons que le modèle devienne linéaire et que les expressions 
correspondant à la sortie soient les suivantes : 


aie +K12 U2 


2.3 
V2=Kysu1+K22 U2 ( : 


où K;; signifie le gain statique déterminé entre y; et u;. Ce modèle peut être écrit plus efficacement en 


notation de matrices comme : 


y = [Klu (2.4) 
Pour les systèmes stables, le gain statique du modèle est, 
K = Gh(s)|._ (2.5) 
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Le gain relatif (4:.) entre y. et u est calculé comme suit : 


ô 
(2) = Kis (2.6) 
U/u, 
0y1 Ki2K21 
——] =K Gi — re) 27 
Ga * K11K22 ( ) 
Y2 
À : (2.8) 
11 : 
1 (22) 
11K22 
les autres éléments de gain relatif sont effectivement calculés pour le cas 2 x 2 comme suit : 
A2 = A1 = 1— A1 and A2 = À (2.9) 
Par conséquent, l'ensemble des gains relatifs pour une trame 2 x 2 peut être symétrique et donné 
comme suit : 
À 1-1 
À = ( _ si 2.10 
1-4  Ân Fe 


Cette équation nous permet de constater que la somme algébrique des éléments de la matrice ;; le 
long d’une ligne i ou d’une colonne j est égale à l’unité. 
Pour les systèmes avec plus de variables (de dimension plus grande que 2x2), l'ensemble des gains 


relatifs peut être déterminé à partir de l'expression, 


A=K@H (2TL) 
où « K » est le gain statique , « @ le produit de Hadamard» et spécifie le produit élément par élément. 
Àj = KijHi;. 
«K;; » est l'élément (i,j) de K et « K;; » est l'élément (i,j) de H = (K-1)T, Par conséquent, le RGA 
est calculé comme suit : 
A=K@ (KT (2.12) 


Il est à noter que pour des matrices non-carrées, K est remplacée par son pseudo inverse. 


2.2.3 Interprétation de la matrice RGA 


La matrice RGA peut être utilisée comme un guide conduisant à un choix des paires des 
entrées-sorties où l’effet des interactions entre les différentes variables caractérisant un 
système est minimale. L’examen des valeurs de ces éléments permet de donner une 
information sur le niveau d’interaction entre ces variables. L’interprétation des valeurs 
Gi; 
Ÿ Si les éléments de la diagonale de la RGA (À 


i — j) peut être classée selon les catégories suivantes : 

ist = j) sont proches de 1, alors le niveau 
55 . X à . . . . 

d’interaction dans le système est très faible, dans le cas contraire les interactions sont fortes. 


Ÿ”_ Si les éléments de la RGA sont très grands, cela signifie que le système est très sensible à de 


faibles variations (erreurs) et donc absolument difficile à contrôler. 

Les cas suivants peuvent se présenter [59]: 

1) À; =1: I n’y a aucune couplage entre les variables du système. Ou l'interaction n’est 
pas dans les deux directions . 

2) À; j = 0 : indique qu'il n'y a aucun effet de l'entrée u; sur la sortie y; . 

3) 0.5 < À;; < 1 : Malgré la présence de l’interaction entre les boucles de régulation correspondantes, 
cette sélection est préférable car elle réduit au minimum les interactions. 

4) À;j = 0.5 : Il y a un niveau important d'interaction. Les autres boucles de régulation ont le 
même effet de l’entée u; sur la sortie y; . 

5) À;j>1 : Il existe un fort couplage entre les boucles de régulation, il doit donc être affaibli. 
L'instabilité peut se produire si À;; est trop élevé. il n’est pas recommandé de contrôler y; avec 
u;. 

6) ) À;; <0 : il y a de fortes interactions, la réponse de la boucle correspondante peut 
changer de sens de variation (système à réponse inverse), si les autres boucles sont fermées. 
En plus la boucle elle-même peut être instable ou le système global devient instable si jamais 
la boucle considérée s’ouvre, d’où le couple correspondant ne doit pas être choisi dans la 


configuration de commande. 
2.2.3.1 Exemplel 


(2.13) 


= Fe | _ LÉ ed 


K21 K22 2 
La matrice RGA est déterminée par l'équation (2.12). 


_ [2.29 —1.29 
— |lEt29 2:29 


L'appariement recommandé est y. avec u, et y, avec u1. 


À (2.14) 


2.2.3.2 Exemple2 


K = ee ci (2.15) 


0.64 0.36 
0.36 0.64 


L'appariement recommandé est y; avec u: et y, avec u2. 


A=| (2.16) 


2.2.3.3 Exemple3 
La matrice de gain à l’état stationnaire pour un système 3x3 est donnée par : 


—2 15 1 
K=1]15 1 —2 


1 —2 135 


(2.17) 


La matrice H = (K-1)T peut être trouvée comme suit: 


—2 15 1 
He(k Tes. À 27 (2.18) 
1 —2 15 
la matrice de gain relatif est donnée par : 
0.7442 1.1860  —0.9302 
A=1-0.9302 0.7442 1.1860 (2.19) 
1.1860 —0.9302  0.7442 


L'appariement recommandé est y. avec u;,y, avec u, et y; avec u,. parmi les six configurations 


possibles. 


2.2.4 Stabilité de la structure de commande sélectionnée 


Même si la méthode RGA est habituellement suffisante dans la plupart des cas, elle ne 
considère pas la stabilité de la structure de contrôle qui en résulte. Aïnsi, pour analyser cette 
stabilité, l’Indice de Niederlinski (ND) est un moyen souvent utilisé [3] Une valeur négative de 
cette indice indique une instabilité dans la boucle de commande [3]. Son expression est 
donnée par l’équation suivante : 

det[G(0)] 
7 I gi CO) < 


L'objectif des deux analyses (RGA et NT) est de déterminer quantitativement la combinaison 


(2.20) 


qui permet de donner un couplage optimal entre les variables d’un processus donné. Ce 
couplage optimal peut être obtenu à partir des deux règles suivantes : 

Règle 1 : L’élément positif de la matrice RGA les plus proches de l’unité devrait avoir 
l’appariement des variables manipulées et contrôlées qui leur correspondent. 

Lorsque l’indice NI est grand, cela implique que le système est moins découplé. 

Règle 2 : Si l’indice NI est négatif, tout appariement entres les entrées-sorties du système est 


inacceptable. 


2.3 Découplage des systèmes multivariables (2x2) 


Lorsque les interactions analysées par l'approche RGA indiquent qu'il existe des 
interactions fortes entre les boucles de régulation, dans ce cas la commande décentralisée 
n'apporte pas de réponses appropriées. Pour surmonter ce problème dans le système de 
commande MIMO), il est souvent préférable d'utiliser une autre procédure qui d'éliminer 
complètement les effets de l'interaction de la boucle est appelée : le découplage [60-61 -62- 
63-64]. Ce qui nous permet le passage d’une synthèse multivariable à un ensemble de 


synthèses monovariables pour déterminer séparément les paramètres du contrôleur pour les 


boucles découplées [59]. Le présent schéma (figure 2.3) montre la structure générale de la 


conception du découplage pour un système multivariable ayant deux entrées et deux sorties. 


Figure 2-3:Système multivariable 2x2 avec découplage 


La figure. 2.3 montre un système de contrôle de découplage pour un procédé multivariable 
2x2. Les variables r,; et r, sont les points de consigne, c; et c, sont les sorties du contrôleur, 
u, et u, sont les variables manipulées et y. et y, sont les sorties du processus. 

En insérant la matrice de transfert de découplage (s) , entre le système MIMO (2X2) et le 
contrôleur, tel que D(s) en série avec G(s), il en produit une matrice de transfert Q(s) 
diagonale, Ou en d'autres termes, il en résulte deux boucles découplées. 

Les éléments diagonaux résultants de Q(s) = G(s)D(s) = diag(q1(s), g22(s))sont : 


CHOSE | 
0  g,,(s) 


où: Les matrices de transfert G(s) et D(s) sont définies par les équations (2.22) et (2.23), 


Q(s) = G(sS)D(s) = | (2.21) 


respectivement : 
. G11(S) G12(5) 
G(s) = RE (2.22) 
et 
D(s) = D:11(S) D:2(s) (2.23) 


D21(S) D22(5) 
Les systèmes découplés (q11(s),q22(s)) peuvent être contrôlés à l'aide des contrôleurs de 
boucles indépendants (ou contrôleurs diagonaux) appelés : contrôleurs décentralisés. 


La matrice de transfert C(s) du contrôleur est diagonale et donnée comme suit : 


Ci1(s) 0 | 
0 C22(s) 


Trois techniques de découplage d’un système multivariable sont dérivées de ce principe 
(Shinskey, 1996 [65]): 


- Découpleur idéal, 


Gt) = | (2.24) 


- Découpleur simplifié, 
- Découpleur inversé 
où le calcul de la matrice de transfert D(s) joue un rôle primordial dans le découplage d’un 


procédé multivariable. 


2.3.1 Découplage idéal 


La substitution des équations (2.22 et 2.23) dans l'équation (2.21) conduit à : 


adj(&(s)) je 0 | 
det(G(s)) 0 q22(S) 


où: adj(G(s)) et det(G(s)) désignent respectivement l’adjoint et le déterminant de la matrice 


D(s) = G(s) * Q(s) = (2:25) 


G(s)donnés par : 
G22(S) —G12(S) 


CE —G21(S)  G11(5) (2:26) 

det(G(s)) = G11(S)G22 (s) = G12 (S)G21(5) (2:27) 
La substitution des équations (2.26 et 2.27) dans l'équation (2.25) donne : 

_ _1 …: 1 G22(5)Q11(S)  —G:12(S)q22(S) 

oo, G11(5)G22(S) — G12(5)G21(8) L-G21(s)quts)  G11(S)Q22(s) | sis 

On peut réécrit l’équation (2.28) comme suit : 
G22(S)411(s) —G12(5)q22(s) 
_ [Gu(s)G22(s) — G12(5)G21(5)  G11(S)G22(s) — G12(5)G21(S) 
PE —G21(5)q11(S) G11(5)q22(s) en 


G1(s)G22(5) — G12(5)G21(S) G11(5)G22(5) — G12(S)G21(s) 


Les éléments G,,(s), G:2(s), G21(s) et G2(s) de l'équation (2.28), qui représentent les fonctions 
de transfert du procédé, sont supposés connus. Les seuls éléments inconnus sont q:,(s) et 
q22 (5). Ils représentent la dynamique désirée du système découplé. 

Dans ce type de découplage, puisque Q(s) est une matrice de transfert diagonale, le choix 
logique pour les composantes de Q(s) est Q.,,(s) = Gas) et Q,(s) = G2(s). Avec ce choix, 
appelé découplage idéal, conduit d’obtenir une matrice de transfert diagonale simple où les 
boucles entrées-sorties sont totalement indépendantes. Les éléments des contrôleurs 
diagonaux Ci1(s)et C22(s) sont réglés indépendamment et respectivement en fonction des 
fonctions de transfert Q,.(s) et Q(s) du système découplé. Un avantage d'une conception du 
découplage idéal est que les contrôleurs C(s) peuvent être réglés facilement et 
indépendamment du fait que l’une ou l’autre boucle soit en mode automatique ou manuel. 


Cependant, cette technique est conduite souvent à des expressions complexes de D(s), qui 
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peuvent être difficiles à réaliser. 


2.3.2 Découplage simplifié 


Une seconde technique de découplage, appelée « découplage simplifié » donnée par 
Luyben [66] qui facilite sa réalisation, est largement utilisée dans la littérature. Elle consiste à 


sélectionner le découpleur comme suit : 


: 1 Du(s) 
ns=|, 4» . | (2.30) 


La figure 2.4 montre l'utilisation de la structure de découplage simplifié d'un système 


multivariable: 


boucle 1 


Figure 2-4:Schéma fonctionnel d’un système multivariable 2x2 découplée par Découplage 


simplifié 


A partir des équations (2.21) et (2.31), les éléments des matrices de transfert D(s) et G(s) 


sont calculés comme suit : 


of ELA, MÉEO 0] en 
L'équation (2.31) permet d'obtenir l'ensemble d'équations suivantes : 
G11(5)+G12(5)D21(s) = Qui(s) (2.32) 
G11(S)D:12(s) + G12(s) = 0 (2.33) 
G21(S)+G22(s)D21(s) = 0 (2.34) 
G21(5)D12(S) + G22(S) = Qz2(s) (2.35) 


On peut dériver les expressions D:32(s) et D21(s), des équations (2.33) et (2.34) 


respectivement comme suit : 


1 
: 1 D:2(5)] _ G11(S) 
ee ES 1 | E . G21(5) 1 16) 
G22(s) 
Alors la matrice de transfert du système découplé est : 
G G 
G11(S) — ere m 2 0 
En G21(S)G12(S) Lens 
0 G22(s) — ta 


Ce choix facilite la réalisation du découpleur simplifié, mais la matrice de transfert du 
système découplé Q(s) obtenue est complexe puisque ses éléments sont la somme de 
fonctions de transfert. Le réglage du contrôleur peut donc être difficile. Il est alors souvent 
suggéré d’approximer les éléments diagonaux par des fonctions de transfert plus simple pour 


faciliter le réglage des contrôleurs indépendants. 


b) 22 


Ya 


Gz1 (s)G42 (s) 


q22(5) = Gz2(s)— Gas) - 


Figure 2-5:Systèmes TITO découplés 


Le processus donné G{s) est découplé à l'aide de la matrice de découplage simplifiée Pour 
concevoir la matrice de découplage, deux cas sont considérés ici. 

Cas 1: Si Les éléments hors diagonale de G(s) n'ont pas de pôles RHP (right-half-plane 
(RHP) poles) et les éléments diagonaux de G{s) n'ont pas de zéros RHP (right-half-plane 
(RHP) zeros) , La matrice de découplage simplifié est: 


L v1(s) Di2(s)v2(s) 
O= Que nt Fe 
Où: 

fa T21 2 T22 
1 (s) = Le. Pre 


1 Ti2 2 Tai 
= 2.39 
va (s)  . Ti2 < Ty ( ) 
gi2(S) 4 
D (s) = — ———— & (Ti2-Tan)S 
" gs) 


g21(S) 
g22(s) 
Cas 2 : Si Les éléments diagonaux n'ont pas de pôles RHP et les éléments hors diagonale n'ont 


e7(t21-T22)S 


D31(S) = — 


pas de zéros RHP de G{s). La matrice de découplage simplifié dans ce cas est: 
Du(s)vs(s) v3(s) 


: _ 2.40 
(s) | va (s) D2(s)w4(s) | 
Où: 

L T22 2 T21 
U3 (s) _ Le T22 < T21 

1 Ti1 2 Ta2 
_—. . nn (2.41) 
92) {4 
Di1(s) = — e”(f22- 720$ 
- 921 (s) 


PACS DR SRE 
D (s) = — ——— € (T11—T12)S 
ui g12(s) 


2.3.3 Découplage inversé 


Shinskey [65] a proposé une technique de découplage intéressante appelée découplage 
inverse. Cette technique consiste à modifier la structure de découplage idéal de la figure 1. 
D'après cette figure et en tenant compte l'équation (2.29), lorsque @Q,,(s) = G1(s) et @22(s) = 
G22(s) , la relation entre les sorties U; du contrôleurs et les sorties U;du réseau de découplage 


est donnée par : 


G22(S)q11(S) —G12(5)q22(S) 
:| — G11(5)G22(5) — G12(5)G21(5) G11(5)G22(5) — G12(5)G21(5) 5 (2.42) 
U; —G21(5)q11(s) G11(S)422(s) U} ‘ 
G11(S)G22(5) — G12(5)G21(5) G11(5)G22(S) mn G12(5)G21 (5) 
Finalement, les éléments de la matrice de transfert de découplage inversé sont 
G12(S) 
1 = Ui — U; 2.43 
1 1 2 Gus) ( ) 
= . G21(S) 
U, = U — U; Gt) (2.44) 


la structure de découplage inversé et indiqué par les équations (2.43) et (2.44) peut être 


représenté comme l’illustre la figure 2.6. 


Figure 2-6:Un schéma fonctionnel d’un système multivariable 2x2 découplée par 


Découplage inversé. 


La structure de découplage inversé évite les inconvénients de découplage conventionnel. 

La structure de découpleur inversé est simple, le système découplé obtenu est simple et la 
synthèse du contrôleur sera aussi simple. 

Malgré les avantages de la technique du découplage inversé, le découpleur ne soit 

réalisable que lorsque ses éléments soient propres et stables. Selon que le système à découpler 


présente ou pas des temps de retard et/ou à non-minimum de phase, le calcul direct des termes 


qe et ee peut conduire à un découpleur avec des éléments ayant des pôles instables. 
11 22 


Ainsi, pour faire face à ce genre de situation, des mesures essentielles de stabilité et de 


réalisation du découpleur doivent donc être prises. 


2.3.4 Découplage statique: 


Il n'est pas toujours possible d'obtenir un découpleur stable. Dans ce cas, il existe d'autres 
manières alternatives de résoudre ce problème de stabilité. Parmi eux, il y a le découplage 


statique qui est donné par les équations suivantes [67] : 


Da2(s) = — lim _ (2.45) 
L 4 G21(S) 
D,,(s) = —-lim Ge (s) (2.46) 


Parmi les avantages du découpleur statique, on peut citer entre autres : la possibilité de sa 
réalisation et la facilité de sa mise en œuvre pratique. 
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2.4 Découplage des systèmes multivariables (3x3) 


On considère la matrice de fonction transfert suivante d’un système ayant trois entrées et trois sorties 
y y 


[68-69-58] : 


G11(S) G32(5) cn) 
(2.47) 


so fente G22(S) G23(s) 
G31(S) G32(5) /Gs3(s) 


Le cofacteur de F(s) = (adj Gp (s)) et la matrice de découplage correspondante sont obtenus comme : 


Fats) F2(s) F3) 
F1(s) Es) P3(s) 
F1(s) F32(s) F33(s) 
(Gas (s)G22(S) — Gaz (s)G23(5)) —(G33(5)G21 (8) " G31(5)G23(9)) —(G31(5)G22 (8) En G32(5)G21(S)) 
—(Ga3(S)G:2 (s) — Ga2(S)Gia(s)) (G33(5)G1(s) — G31(5)G:3 (s)) —(G32(5)G:1 (5) — Ga1(S)Gs2(s)) (2.48) 
—(G11(5)G1(s) : G1(5)G(s)) —(G11(5)G(s) — G11(5)G1(s)) (G11(5)G1 (8) mn Gu(s)Ga(s)) 
F1(s) F31() 


F(s) = 


1 ue 
a. F22(5) sn 1 D;2(s) Di3(s) 

D(S)simptité = | 5 rt CO 0) (2.49) 
F1(s) F3:(5) Das (s) D: (s) 1 


F3(s) P3(s) 
F1(s) F:2(5) 


Les éléments du découpleur sont : 
G33(S)G12(S) — G32(S)G13(S) 


1 


G22(5)G13(S) — G23(S)G:12(s) 


Dis) = G33(5)G11(S) — G31(5)G13(5)” PGI G22(5)G11(5) — Gz21(5)G12(S) 

node G33(5)G13(5) — G31(5)G23(5) D G23(S)G11(s) — G21(s)G3(s) 
G33(5)G22(5) — G32(5)G3(5) G32(5)G11(S) — G21(5)G12(s) 
G31(5)G22(5) — G32(5)G21(s) __ G32(5)G11(s) — G31(s)G12(s) 

Die, 06 


” Gas 65) Gals) | OMC EMEA) 


2.5 Conception des correcteurs conventionnelle (PI/PID) des systèmes 
découplés 

La conception du contrôle des systèmes multivariable est généralement divisée en deux 
types : la commande centralisé et le la commande décentralisé [70-71-72-73-74]. La 
commande décentralisée est la plus utilisé dans l’industrie, vu sa simplicité de mise en œuvre, 
sa robustesse et son flexibilité. Ainsi les performances peuvent être obtenues en se 
concentrant sur les performances des boucles individuelles. 
Parmi ces avantages, on peut citer : 
1. la simplicité de mise en œuvre pratique et d’entretien, 
2. le maintien de la stabilité dans le cas où l’une des boucles se trouve coupée (défaillance 


d’un capteur ou d’un actionneur), 


3. la non-propagation de la perturbation, agissant sur une sortie, dans le système, 

4, la possibilité d’application des techniques de commande utilisées en monovariable et les 
performances qu’on peut atteindre en s’intéressant aux performances de chaque boucle 
indépendamment de l’autre, 

La conception des contrôleurs décentralisés nécessite également le réglage des paramètres de 
ces contrôleurs. Ce réglage se fait d’une façon indépendante d’une boucle de régulation à 
l’autre, permettant de contrôlées les entrées et sorties de manière indépendante. Si les 
systèmes découplé sont d’ordre élevé dans ce cas, il est difficile de concevoir directement les 
contrôleurs pour les éléments diagonaux de matrice (s). Et donc, plusieurs techniques 
différentes ont été proposées pour approximer les éléments diagonaux ((Q:1(s),Q22(s)) par 
des modèles du premier ou deuxième ordre avec retard (FOPDT/SOPDT) [75-76-71]. 


où: FOPTD (en anglais : First Ordre Plus Time Delay), SOPTD (en anglais : Second Ordre Plus Time 
Delay). 


2.5.1 Modélisation des systèmes découplés par une fonction de transfert d'ordre entier 


2.5.1.1 Méthodes de Wang 


Dans cette section, la fonction de transfert exprimée dans le domaine de Laplace de tels modèles 
Gu(s)i de chaque élément Q(s);; de Q(s ) s’écrit [77-78]: 


Ki 


——— Le-Lis il = 1,2. 2.50 
Ts + 1° ; É ( ) 


Gu(s)ii = 


où: K;, est le gain du modèle ; L;;est le temps de retard et leT;; est la constante de temps 
Trois paramètres inconnus (K;, L;; et T;;) dans l'équation. (2.50), qui doit être spécifiée, afin 
d'obtenir le modèle FOPDT de Q(s);;. Dans cette section, un modèle FOPDT a été obtenu sur la base 
de la réponse en fréquence appropriée en deux points, w = Det & = wc;; où wc;; est la fréquence 
d'intersection de phase [78]. 
Gn(O)i = Q; (0) 


1Gu Waii)l = [Q. GWi)| (2.51) 


LGmOWoii) = 20, OWaii) 


En conséquence, les paramètres du modèle FOPDT peuvent être calculés en utilisant 


Ki = Qi (0) 


KF — Qi GWe)l2 
Qui OGwei)l we; (2.52) 


_T+ tan !(-w,iTi) 
WoiiTü 


ii 


2.5.2 Méthodes de synthèse du régulateur PI/PID décentralisés 


La stratégie de réglage des correcteurs monovariables (ou synthèse monoboucle de 
correcteurs PI/PID) est importante dans le cas d’une stratégie décentralisée. Divers types de 
méthodes classiques de synthèse monoboucle permettent le réglage de correcteurs PI/PID. Les 
principales sont les méthodes heuristiques [79-80], les méthodes analytiques [81-82], les 
méthodes d’optimisation [83-84] et les méthodes de modelage [85-86]. 

La plupart des méthodes de synthèse multiboucles s’appuient sur des méthodes de synthèse 
monoboucle de correcteurs PI/PID. 

On distingue quatre types de synthèse multiboucle : 

- Les méthodes indépendantes [87-88]. 

- Les méthodes séquentielles ou itératives [89-90]. 

- Les méthodes dites simultanées [91-92]. 


- Les méthodes de detuning [93-941]. 


Cependant, les performances résultantes peuvent être médiocres car les méthodes 
manquent d'informations détaillées sur la dynamique du contrôleur dans d'autres boucles. 
Toutes les méthodes de synthèse multiboucles disponibles dans la littérature ont tenté de 
surmonter l'interaction entre les boucles. Le réglage des paramètres du contrôleur d'une 
boucle affecte les performances des autres et conduit parfois à une instabilité de l'ensemble du 
système. Pour des performances stables, de nombreux contrôleurs décentralisés industriels 
sont réglés de manière lâche, ce qui entraîne un mauvais fonctionnement et des coûts 
énergétiques plus élevés. Un autre moyen possible de surmonter cette difficulté et d'utiliser les 
méthodes de conception SISO est de construire des boucles individuelles équivalentes à l'aide 
d'un découpleur. Les contrôleurs PI/PID décentralisés avec découpleur sont le choix le plus 
courant pour les systèmes (2x2). 


Dans cette approche, plusieurs méthodes de réglage pour les contrôleurs PI/PID 


décentralisés pour les systéme multivariables basées sur le découpleur ont été proposées, 


parmi lesquelles: 


2.5.2.1 Méthode NDT 
Les paramètres du régulateur PI par la méthode NDT ((Non-Dimensional Tuning method [95 NDT]) 


sont déterminés en fonction des paramètres du modèle obtenu dans l'équation (2.50), comme indiqué 


Ty 1 
er = 2.59 
pii Kii 21; Ba 14 ( ) 
Ti ( 1 Li; si) 
ini 60 2.54 
mn V7 n 14 ee) 


2.5.2.2 Méthode SIMC 


Les paramètres du régulateur PI par la méthode SIMC ( Simplified Internal Model Control [96 


SIMC)) sont: 
À 5%, 
Ru = (2.55) 
PE UK Toi + Li 
Ti = min(Ti, (Toi + Li)) (2.56) 


Pour plus de détails sur les méthodes de conception NDT et SIMC (voir dans [95-96]). 


2.5.2.3 Méthode GPMs 


Les paramètres du régulateur PI par la methode de synthèse GPMSs (gain and phase margin 


specifications [77]) sont données par : 


WpiT ii 
Ki Ta (2.57) 
1 
4w2;Li 1 
Ti = (2W4 — + — 2.58 
Ii | Wii T +) ( ) 
et: 
1 
Asnii Omii + 3 T(Amii — 1) 2 59 
PET (Ami — 1) Li en 
OÙ: Ki , Ta, Lu les Paramètres du modèle obtenu dans l'équation (2.50), 
Amii Et Oynii les valeurs désirées de marges de gain et de phase respectivement. 
2.5.2.4 Commande par modèle interne (Internal Model Control). 
Les paramètres du régulateur PI sont [97]: 
2Ti + Li 
Koss = ————"— (2.60) 
FF 2Ku(Lu + B) 
1 
Ki ==> > (2.61) 


Ki (Li + B) 


Le paramètre £ devrait satisfaire : 8 > 0.2 Ti;et B > 0.25L;;. 


2.6 Application : 


On considère la matrice de fonction transfert d'un modèle de la colonne de distillation présentée 


par Wood et Berry (1973[98]) ayant deux entrées et deux sorties (2x2) est écrit sous la forme : 


12.8e 15 —18.9e 3 
16.75 + 1 21s+1 
6.6e 7  —19.4e 55 
10.9s+1 14.45 +1 
Les éléments de la matrice RGA et la valeur de l'indice NI sont : 


2.0094 — 1.0094 
RGA(G(O)) = —1.0094 2.0094 l 


G(s) = 


NI=0.4977 


Nous observons à partir de l'analyse de la matrice RGA et de la valeur de NI que le meilleur 


appariement recommandé est (u1-y1) ef (U2-V2). 


boucle 1 


.477(16.75+1) _,. 
215 +1 
0.34(14.45+1) _,. 
— 
10.95 +1 


16.75 + 


Figure 2-7:Commande décentralisé d'un système découplé 


La matrice de fonction de transfert du découplage simplifié déterminé à l'aide de l'Eq. (2.30) est: 


1.477 (16.75 + 1)e”?2s 
215 +1 


4s 1 


1 
0.34(14.4s +1) _ 
10.9s + 1 
Les éléments diagonaux résultants de Q(s) = F(s)D(s) = diag(q11(s), q22(s))sont : 


Di -28e" | 643(1445+ 1e 7 
MS 67s+1 (10.95+1)(215+1) 


D;(s) = 


— 19467, 9.745(16.75 + 1)e*° 
14.45 +1 (10.95+1)(21s +1) 
Ainsi, en utilisant l'équation (2.66) , Les modèles FOPDT réduits de h:1(s),h22(s) sont les suivants : 


À D- 7 
M1 FOPDT ST 1s+1 


Q22(S) = 


—1.0655s 


—9.655 


G = —2.157S 
m2 ropr (S) 4.684s+1 € 
Les paramètres de régulateur PI par la méthode NDT: 
0.074 
0.41 + nan 
C(S)nor-p1 = 0.024 
(à) —0.12 — — 
Les paramètres de régulateur PI par la méthode BLT: 
0.045 
0.375 + —. (à) 
C(S)gir-P1(S) = 0.003 
(à) —0.075 — — 
s 
Les paramètres des régulateurs PI et PID par la méthode GPMSs sont donnés respectivement : 
0.0881 
0.4867 + na 
GcMaghade,P-P1(S) = 0.0304 
(à) —0.1567 — 
s 
et 
0.9733 + 0.0887 + 2.68875s 
(à) 
_ 9.5225 +1 
GcMaghade,P-PID (s) — 0.0304 
—0.3134 — = — 0.8070s 
(à) 
5.1499s +1 
Les paramètres de régulateur PI par la méthode SIMC: 
0.0737 
0.3988 + EE (à) 
C(S)simc-p1 = 0.0240 
0 —0.1125 — 
Les paramètres de régulateur PI par la méthode IMC: 
0.0177 
0.3054 + ic (à) 
C(S)imc-Pi(s) = 0.0054 
(à) —0.0863 — 


Le temps de simulation totale estimée est 200 secondes. Un échelon unitaire est introduit aux instants 
t-0 et t-100 pour les références des boucles respectivement (u:-y.) et (u2-y2). Les résultats de la 


simulation sont présentés dans la Figure. 2.8. 


—— GPMs-PID 
ze l |—— GPMs-PI 7 
0.4 — NDT-PI | 
_ — SIMC-PI | 
—— IMC-PI 
ES RS RS RS RS RE — BLT-PL 1 
-0.2 ; | | | 
[el 20 40 60 80 100 120 140 160 180 200 
Temps (sec) 
(a) 
1.6 T T T T ï T T 1 
14 L ——— GPMs-PID | 
sl |—— GPMs-PI | 
— NDT-PI 
LF —— SIMC-PI 
__ [—— [MC-PI 
>=" ——— BLT-PI 


-0.2 1 i i i i i ñ i 
[e) 20 40 60 80 100 120 140 160 180 200 


Temps (sec) 


(b) 


Figure 2-8:Réponses indicielles du système contrôlé (Système MIMO 2 x 2) avec les 
contrôleurs classiques 


La figure 2.9 montre les signaux de commande appliqués au système étudié 


—— GPMs-PID 
—— GPMs-PI 
—— NDT-PI : 
— $SIMC-PI 
—— IMC-PI 
— BLT-PI 


-0. 1 L Ï 
O 20 40 60 80 100 120 140 160 180 200 
Temps (sec) 
(a) 
0.1 T T T T T T T T T 
0.05 


-0.05 — GPMS-PID 
.. — GPMs-PI 
-0.1 +  NDT-PI 
—— SIMC-PI 
Sos | —— IMC-PI 
_ BLT-PI 
-02 + 
-0.25 
O 20 40 60 80 100 120 140 160 180 200 
Temps (sec) 
(b) 


Figure 2-9:Signaux de commande classiques appliqués au système multivariable. 
2.7 Conclusion: 


Dans ce chapitre, Une étude générale des systèmes multivariable a été présentée, en 
commençant par l'identification des techniques permettant de surmonter les problèmes 
d'interférence entre les entrées et les sorties dans les boucles de régulation, puis nous avons 
présenté les trois méthodes de découplage les plus répandues dans la littérature ( découplage 
simplifiée, découplage inversé et idéal) des systèmes multivariable de type 2X2 et 3x3. 
Plusieurs méthodes de synthèse des contrôleurs PI décentralisés classique combinés avec 
découplage simplifiée ont été utilisés commander la colonne de distillation comme exemple 


d’application. 


Chapitre 3 : Synthèse optimale des contrôleurs fractionnaires 
pour les systèmes multivariables 


3.1 Introduction: 


Le but de ce chapitre est de trouver des techniques plus efficaces et optimales pour 
déterminer les paramètres du contrôleur d'ordre fractionnaire afin d’imposer les performances 
désirées aux système smultivariables carrée. Nous avons utilisé les algorithmes 
d'optimisation d'essaims particulaires (PSO: Particle Swarm Optimization) pour déterminer 
les paramètres optimaux en minimisant une fonctionnelle du coût en fonction de l’erreur de 


type ISE, IAE et ITAE. 


3.2 Optimisation par essaim particulaire (PSO) : 


L’optimisation par essaim particulaire (PSO), est un algorithme évolutionnaire qui utilise une 
population de solutions candidates pour développer une solution optimale au problème. Cet algorithme 
a été proposé par Russel Eberhart (ingénieur en électricité) et James Kennedy (socio-psychologue) en 
1995 [99] Il s’inspire à l’origine du monde du vivant, plus précisément du comportement social des 
animaux évoluant en essaim, tels que les bancs de poissons et les vols groupés d’oiseaux. 
Dans l'algorithme PSO de base, chaque particule : a deux équations de base qui sont la mise à jour 
de la vitesse et la mise à jour de la position comme présenté ci-dessous [99-100] : 
e La vitesse : 
Vo = WG — 1) + cri (Pipese( — 1) — Xi — 1)) + cor (Gijpese(t — 1) — Xi —)) (3.1) 
avec :jEe{1..D} ,ie{1...P} 
e La position: 
Xi, CE) = Xij(t1 — 1) + V;;()j € {1...D} (3.2) 
Où: 
Vij (ti — 1) : est la vitesse actuelle de la particule. 
Vi; (E) : est la vitesse mise à jour de la particule. 
Xij(t1 — 1) : est la position actuelle de la particule, 
ZX; (€): est la position mise à jour de la particule. 
rl etr2 : sont des nombres aléatoires compris entre {0,1}. 
c1 : est le coefficient cognitif (généralement 2,0) et c2 : est le coefficient social (typiquement 1,8). 
w : est le coefficient de masse d'inertie (généralement 0,7). 


Le coefficient d’inertie peut être variable selon [101-102] : 


W = Wmux — (as tr) x k (3.3) 


Kkmax 


OÙ: Kmax ; À Sont respectivement le nombre maximum des itérations et le nombre d’itération courante. 


Wmax » WminSOnt respectivement les coefficients minimum et maximum d’inertie. 


3.2.1 Principe de fonctionnement de PSO 


La position de la particule est mise à jour en utilisant sa position et sa vitesse actuelles. La mise à 
jour de la vitesse est basée sur la meilleure position de particule mémorisée (meilleure position locale) 
et la meilleure position de ses voisins (meilleure position globale), de sorte que les ajustements de 
position de chaque particule seront basés sur l'apprentissage historique de sa trajectoire et des 
meilleures positions voisines. 

Soit /{x) la fonction objective à optimiser (fitness) et 7 le nombre de particules. Les étapes essentielles 
de l'optimisation par PSO sont présentées par l'algorithme suivant : 

Etape 1: initialiser une population de particules et de vitesses d’une façon aléatoire, uniformément 
reparties à l'intérieur de l'espace de recherche D prédéterminé, et fixer le critère d'arrêt. 

Etape 2 : évaluer la fonction objective, pour chaque particule. 

Etape 3 : mise à jour de la vitesse : à chaque itération, les vitesses de toutes les particules sont mises à 
jour selon l’équation (3.1). Si la nouvelle valeur de vitesse de particule excède les limites 
prédéterminées, appliquer la stratégie de confinement de la vitesse. 

Etape 4 : mise à jour de la position : à chaque itération, les positions de toutes les particules sontmises 
à jour selon l’équation (3.2). 

Etape 5 : mise à jour du coefficient d’inertie : à chaque itération, le coefficient d’inertie est mis à jour 
selon une loi décroissant de l’équation (3.3). 

Etape 6 : mise à jour des meilleures particules : à chaque itération, mettre à jour la meilleure position 
atteinte personnellement (Pest) et la meilleure position atteinte par l’essaim (Gpest) lorsque la 
condition est remplie. 

Etape 7 : vérifier le critère d'arrêt : répéter les étapes 3 à 6 jusqu'à ce que le critère d'arrêt (nombre 
maximum d'itérations atteint) soit satisfait. 


L’organigramme de la norme PSO est présenté par la Figure 3-1. 


Initialisation aléatoire de position 


et de vitesse des particules 


Mettre à jour la vitesse des particules 


Confiner La vitesse des particules 


Mettre à jour La position des particules 


Mettre à jour Le coefficient d'inertie 


Critère 
d'arrêt ? 


Non 


- » 


Meilleurs paramètres 


Figure 3-1:Organigramme de l’algorithme de PSO utilisé. 


3.2.2 Optimisation des paramètres optimaux du contrôleur fractionnaire 


Le principe générale de notre recherche est de trouver des valeurs optimaux pour les paramètres du 
contrôleur PID d'ordre fractionnaire par l'algorithmes d’optimisation PSO, un problème d'optimisation 
est formulé en utilisantles formules d’erreur suivant : ISE, IAE et ITAE, La Figure 3-2 montre la 


structure générale de la méthode proposée. 


Figure 3-2:Principe de l’optimisation par algorithme PSO. 


Les cinq paramètres K,,, K;, Ka, À,u sont recherchés et dans le domaine de stabilité et ceci en 


faisant le choix selon un critère d’optimisation J définie par [103] : 
hse(K, Ki, Ka, À, 4) = | e(t)? dt 
0 


Je (Kps Ki, Kay À, Li) | le(t)l 
0 


= (34) 


Jiraë (K, Ki Ka, à li) =) tle(t)|dt 
0 


Jirse (Ko, Ki, Ka, À,U) = :| e(t)? dt 
0 


3.3 Procédures de conception pour les systèmes multivariables 


Dans notre travail, l’algorithme PSO est utilisé pour déterminer les paramètres optimaux des 


correcteurs PI d'ordre fractionnaire décentralisé. 


3.3.1 Premier cas (un système à deux entrées et deux sorties (2x2)) 


Pour système de deux entrées et deux sorties (2x2), illustré à la Figure 3-3. 


Paramètres initial 


Paramètres Optimisés 
minJ(K,.K:7%,.7: À À) 


découplage (2x2) Procédé (2x2) 


ue) (OS 620) 
Gz1(s) Gz2(s) 


Dz1(s) Dz2(s) 


Figure 3-3:Principe d’ajustement des paramètres du régulateur PI d'ordre fractionnaire par 


l’algorithme PSO. 


Les systèmes découplés sont contrôlés par deux contrôleurs PI fractionnaires PI et PI/2 
(Figure 3.3). Les fonctions de transfert des deux contrôleurs sont données par les 


équations (3.5 et 3.6). La réalisation diffusive de chaque contrôleur peut être exprimée comme suit. 


FO PILE Î LE dE (3.5) 
FO — PI, = Î ET (3.6) 


Nos objectifs sont de définir les six paramètres optimaux (K51, Ti , Kh2, Ti2 1,42) des contrôleurs PI 
fractionnaires décentralisés dans la région de stabilité et ceci en faisant le choix selon un critère 


50 


d'optimisation Jmin définie par: 
haerse rar (Kp1e Tia > Kp2, Ti2 A1) = 1 f(e(t))dt (3.7) 
0 
Cette fonctionnelle est constitué des indices de performances de type ISE,IAE et ITAE. 


Ÿ__ Critère Integral of Square Error (ISE) : représente l'intégrale de l’erreur quadratique. 


Le critère de performance est présenté par la formule suivante : 
ISE = Î e(t}? dt (3.8) 
0 


Ÿ__ Critère Integral of AbsoluteError (TAË) : représente l’intégrale de l’erreur absolue. 


Le critère de performance est présenté par la formule suivante : 


re | lece)| (3.9) 
0 


Ÿ__ Critère Integral Time multiplied by AbsoluteError (ITAE) :représente l’intégrale de 
l’erreur absolue multipliée par le temps. 


Le critère de performance est présenté par la formule suivante : 


00 


ITAE =) tle(t)|dt (3.10) 
0 


La formulation mathématique de problème d’optimisation est basée sur l’approximation des 
fonctions de transferts irrationnelle des correcteurs fractionnaires par sa représentation 
diffusive implémenté dans Simulink/MATLAB. 
Les étapes pour déterminer les paramètres des contrôleurs PI fractionnaires peuvent être résumées 
comme suit : 
1. Implémenter le système de contrôle de feedback du système contrôlé avec un contrôleur 
PI fractionnaire dans Matlab/Simulink, y compris la réalisation diffusive de l'opérateur 
fractionnaire via le modèle Simulink. 
2. Sélectionnez la fonction PSO appropriée de la boîte à outils d'optimisation Matlab pour 
minimiser la fonction objectif} . 
3. La méthode de réglage non dimensionnel basée sur la marge de stabilité (NDT) [95] est 
utilisée pour déterminer les paramètres initiaux des contrôleurs. 
4. Calculez l'erreur (IAE, ISE et ITAE). 
Plusieurs méthodes d’optimisations ont été proposées dans la littérature, et ont prouvé leur efficacité 
dans différent domaines d’application. Nous nous intéressons dans ce chapitre à une méthode méta 


heuristique, l’optimisation pas essaim particulaire (PSO). 


3.3.2 Deuxième cas (un système à trois entrées et trois sorties (3x3)) 


Pour un système multivariable de trois entrées et trois sorties (3x3) avec simplifies découplage, illustré 


à la Figure 3-4. 


L 


qu 


Figure 3-4:Principe d’ajustement des paramètres d’un régulateur PI d'ordre fractionnaire 


pour un système 3x3 


Les boucles indépendantes découplées sont représentées sur la Figure 3-5. 
Lesboucles décomposées sont contrôlées par trois contrôleurs PI d'ordre fractionnaires PI/1et 
PI2 et PI}: . Les fonctions de transfert des trois contrôleurs sont données par l'équation 


(3.11 )comme suit: 


Il Il Fe 
C(S)11 = k, (pet) co HE [nee ) co T [ns DE 


il i2 É 


La matrice de transfert C OF du contrôleur est diagonale et donnée comme suit : 


C(S)11 0 0 Le 
C(s)ÿ=| O0  C(s)z2 0 |- 0 RUES 0 (3.12) 
0 0 C(s)33 T; 
0 0 (ps) 
l 5 | 


{111421 PSO Alscorithm 


| { » 
Lan iÆS A Æia 7,7, ZA 2 SN) | LE RES 


Figure 3-5:boucles décomposées d’un système ayant trois entrées et trois sorties (3*X3). 


Il s’agit de la recherche des neuf paramètres K,:,Ti:, Kp2, Ti2 K53, Ti3 A1, 2 etAsdes trois 
contrôleurs PI fractionnaires décentralisé dans la région de stabilité et ceci en faisant le choix selon un 


critère d’optimisation J . 


3.4 Application et Simulation 


3.4.1 Application aux systèmes réels (Modèle linéaires de la colonne de distillation) : 


Le modèle mathématique identifié par Vinaya et Arasude colonne de distillation binaire 
d'usine pilote [104] est exprimé par la matrice de transfert suivante (system multivariable ou 


MIMO: de deux entrées et de deux sorties) : 


SOA QE nes 


_[Ti45+1 0645+1 
BOL gages pi8e-00 Fo 


1.235 +1 0.32s + 1 
e Application de découplage simplifié: 


Le calcul des deux termes de découpleurD;, et D;1 pour le système multivariable de la 
colonne de distillation, on trouve : 


1 1.3846(1.14s + 1) 


_ 0.64s + 1 
DS) = |: 888(0.325 + 16-11 i F4 
1.235 +1 
Les éléments diagonaux résultants de H(s) = G(s)D(s) = diag(h:1(s), h2,(s) sont : 
. h11(5) 0 
CD = |  . (3.15) 
où: 
—0.13e-0035  0.3398(0.32s + 1)e- 122 
hu(s) = ER (3.16) 


1.145 +1 (0.645 + 1)(1.23s + 1) 


(= 2186 , 0.47076(1.145 + Ie V2 , 
28 G325+1  (1.235+1)(0.645+1) (3.17) 


3.4.1.1 Synthèse des régulateurs PI décentralisés par les méthodes classique NDT et SIMC: 


Dans ce chapitre, nous sommes appelés à faire une comparaison de les méthodes NDT et SIMC, qui 
sont les méthodes les plus populaire pour le réglage du contrôleur PI décentralisé. 
Les fonctions de transfert des modèles du premier ordre avec retard, qui sont équivalentes aux 
systèmes découplées sont [104]: 

0.21 


= -2s 
Q11 roppr (s) = 1355241 e (3.18) 
et 
—0.29 
— —1.5455s 3.19 
Q11 roppr (S) 0.735s+1 € à ) 
Ÿ”_ Les paramètres du régulateur PI par la méthode NDT sont: 
Ti 1 
Kiki = -+— 3.20 
pui 2L;; 14 ( ) 
Tri ( 1li Li 
F =mi 1+570,98) S21 
nanas TT Fr 
Tableau 3-1:Les paramètres de régulateur PI par la méthode NDT. 
Proportionnel Ki Gain d’intégral K;;; 
Boucle 1 1.9473 1.19 
Boucle 2 -1.0665 -1.116 
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Ÿ”_ Les paramètres du régulateur PI par la technique SIMC sont: 


Ki = Li — (3.22) 
Ki Toi + Li 
Ti = min (Ti, A(Toù + Li)) (3:23) 
Pour plus de détails sur les méthodes de conception NDT et SIMC (voir dans [95-104]). 
Tableau 3-2:Les paramètres de régulateur PI par la technique SIMC. 
Proportionnel Ki; Gain d’intégral K;;; 
Boucle 1 1.6071 1.1904 
Boucle 2 -0.8202 -1.116 


3.4.1.2 Synthèse optimal des régulateurs PI d'ordre fractionnaire décentralisés par la 
méthode proposée 


Les contrôleursPl" ,Pl/2décentralisé est placé en cascade avec la fonction de transfert du 
système multivariable (2X2) avec découplage simplifie donné par le bloc de diagramme dans 


laFigure 3-6. 


| Initial 


l 


parameters Optimization 
(kp1s Kia @1, kp2ki2,@2) 
D(s) . G(s) = diag(h:,(s),hz2(s)) 


E D | e-0.03s 
R1 FO -{PI, 0.12 Yi 


14521 


Y2 
R2 + 4 
nl FO — PI, 


fractional order decoupler x process — decoupledsystem 
controllers 


Figure 3-6:Structure de système de commande par contrôleur PI d'ordre fractionnaire. 


Les paramètres initiaux de l’algorithme PSO utilisés pour optimiser les paramètres des 
contrôleurssont donnés dans le tableau 3.3: 


Tableau 3-3:Paramètres initiaux de l'algorithme 


Paramètres FO — PI, Fractionnaire FO — PI, Fractionnaire 
Taille de la population 20 20 
nombre maximum d'itérations 120 120 
Nombre de particules 6 10 
Limite inférieure et limite —5<K?p1,Kn <5 —5 < Kp2 ,Kp < 5 


supérieure des variables 


0<A <2 0<À <2 
a.est ordre d'intégration aest'ordre d'intégration 
Coefficient cognitif Ci 15 15 
coefficient social C2 2 2 
Coefficient d'inertie W 1 1 


L'objectif de cette partie est d’appliqué la stratégie d’optimisation PSO dont le but de calcul des 
paramètres (Kp1 , Kj1, Kp2 , Kj2, À , A2) Optimaux d’une commande décentralisée. 

Les paramètres optimale du régulateur PI fractionnaire obtenues par l’algorithme de PSO 
avec le critère IAE est donnés dans le Tableau 3.4. 

Tableau 3-4:Les paramètres de régulateur PI/fractionnaire par la technique de synthèse 


optimale propose avec le critère IAE. 


Proportionnel K;; Gain d’intégral K;;; A2 
Boucle 1 1.9278 0.010479 0.13567 
Boucle 2 -0.9786 -0.01038 1.1451 


3.4.1.3 Comparaison des différents régulateurs PI décentralisés 


Les performances du contrôleur PI d’ordre fractionnaire proposé [14] sont comparés avec celles 
obtenus par le contrôleurs PI classique tel que NDT-PI [95-104] et SIMC-PI [104-105], sont discutés à 
l'aide de MATLAB/Simulink. 

Les résultats des performances (IAE, ISE, ITAE, dépassement (%) et Temps de Réponse (sec)) 
obtenus par les deux contrôleurs (d'ordre fractionnaire et contrôleurs PI classique) ont été organisés 
dans le Tableau 3.5. 


Nous comparerons les performances du système de commande avec les deux contrôleurs: le 


contrôleur PI classique (NDT-PI et SIMC-PT) et le contrôleur PI d’ordre fractionnaire. 


Tableau 3-5:Les caractéristiques de la réponse temporelle du système commandé 


Performances Entrée (u)-Sortie Les méthodes de réglage 
(y) Méthode Proposée Méthodes classiques 
FO-PI NDT-PI SIMC-PI 
IAE U1_Y1 5.123 5.946 6.375 
U_ V2 3.55 3.951 4.431 
ISE U1_Y1 4.561 4.797 4.945 
U_ V2 3.725 3.948 4.063 
ITAE U1_Y1 15.06 23.22 28.55 
U_ V2 6.14 8.134 11.44 
Dépassement U1_Y1 6.2 26 27.42 
(% U>_Y2 5.01 23.2 23.95 
Temps de Réponse U1_Y1 8.02 15.12 17.15 
(sec) U>_Y2 7.45 8.11 8.54 


Les figures 3.7 et 3.8 présentent les réponses indicielles des deux boucle (Y,: 1°° boucle 


et Y,:2°% boucle) du système de commande en boucle fermée : 


1.5 T T T T 
CS 
0 
1 LE —. 
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Y1 


=" O-PI Controller 
---NDT-PI Controller 
- - -SIMC-PI Controller 
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Figure 3-7:Réponses indicielles de la première sortie (Y1) du système de commande avec les 


régulateurs PI fractionnaire (FO-PT) et PI classique (NDT et SIMC). 
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Figure 3-8:Réponses indicielles du première sortie (Y2) du système de commande avec les 


régulateurs PI fractionnaire (FO-PT) et PI classique (NDT et SIMC). 


Pour démontrer l'efficacité de la méthode de conception proposée, nous comparons nos résultats avec 


ceux obtenus à l'aide des contrôleurs PI classiques (NDT-PI et SIMC-PT) sont donnés comme suit : 


> 


D’après les résultats obtenus dans les deux figures (3.7 et 3.8) et le Tableau 3.5,0n peut 
déduire que l'algorithme de PSO est très efficace pour optimiser le contrôleur PI fractionnaire 
qui assurer d'excellent comportement fractionnaires du système  multivariable avec des 
performances de contrôle très satisfaisantes. 

Le temps de réponse obtenu par les réponses temporelles du système commandé par le 
contrôleur PI d’ordre fractionnaire est plus petit que celui obtenu par le contrôleur PI 
classique (NDT-PI and SIMC-PI ) avec les trois cas des critères intégraux (ISE, IAE et 
ITAE ). 

On note également que les différents indices de performance (IAE, ISE, ITAE) en tant que 
fonctions objectifs, avec le contrôleur PI d'ordre fractionnaire donnent des résultats très bons 
et supérieurs par rapport aux régulateurs classiques (NDT-PI et SIMC-PD), voir Tableau 3.5. 
Le contrôleur proposé avec le decouplage simplifié est plus sécurisé car le dépassement est 
presque nul, on peut voir ici sa robustesse, les résultats de simulations montrent l’efficacité 
égalemen l'utilisation de découplage simplifié en terme de rejet de perturbation et 
d'atténuation les interactions entre les entrées/sorties du système commandé et les 
performances supérieures du contrôleur PI fractionnaire (FO-PT. 

Malgré la complexité de la structure du contrôleur PI fractionnaire qui a offert en vérité de la 
flexibilité, les résultats obtenus ont été très satisfaites et l’algorithme de PSO utilisé dans 
l’optimisation des paramètres du contrôleur PI fractionnaire est très efficace et donne des 
bonne résultats. En peut donc dire que la performance de contrôleur PI fractionnaire a surpassé 


le contrôleur PI classique ce qu’il lui permet d’ajouter des spécifications additionnelles. 
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3.4.1.4 Tests des changements de consignes 


Dans cette partie, plusieurs tests des changements de consignes ont été faits. Les résultats 
de simulations des réponses temporelles en boucles fermées pour procède MIMO seront 
présentés avec : 


e Dans la première simulation, la valeur désirée de consigne pour la première sortie est 50. et la 


valeur désirée de consigne pour la deuxième sortie est 70. 
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o _l---NDT-PI Controller of Output Y2 el 
SIMC-PI Controller of Output Y1 
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Figure 3-9:Simulations du système contrôlé avec changements de consignes : une consigne 
de 50 est appliquée sur R;, a l’instant t-0 minute. et Une consigne de 70 est appliquée sur R>, 


a l’instant t-0 minute. 


e Dans la deuxième simulation, la valeur désirée de consigne pour la première sortie est 64. et la 


valeur désirée de consigne pour la deuxième sortie est 74. 
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Figure 3-10: Simulations du système contrôlé avec changements de consignes : une consigne 


de 50 est appliquée sur R:, a l’instant t-0 minute. et une consigne de 70 est appliquée sur R>, 


a l'instant t-0 minute. 
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D'après les Figures (3-9 et 3-10), on remarque que nous avons obtenus les meilleures performances du 
système contrôlé avec contrôleur PI d'ordre fractionnaire et n’observe aucune changement 
dépassement, ainsi qu’un temps de réponse court malgré des changements de consignes désirées. On 
peut dire que l'utilisation du contrôleur PI d’ordre fractionnaire décentralisé basé sur l'approche 
diffusive avec découpleur offre des bonnes performances et robustesse. 

On a confirmé dans notre thèse l’efficacité de la commande d’ordre fractionnaire décentralisé FOPI 
appliquée sur le modèle de la colonne de distillation W.B (système MIMO) à travers les changements 


de consignes. 


3.4.2 Application aux systèmes réels (3x3) 
Considérons un processus multivariable à trois entrées et trois sorties donné par la Figure 3.11 
V1 


E 
E 


Figure 3-11:Processus multivariable à trois entrées et trois sorties 


Dans étude de simulation, nous considérons la fonction de transfert d'un processus multivariable 3x3 
donné par [106-107-108] comme suit, 
0.66 —0.61 —0.0049 


e” 2:65 e” 355 e 55 
6.75+1 8.45+1 9,065+1 
1.11 _2.36 _0.01 
G = —6.55s —3s —1.2s 
QE EET-an Le 55s+1° 7.095 +1 
_34.68 on 462 og _ 0.87(11.615+1) 
8.155+1 10.95+1 (3.89s + 1)(18.8s + 1) 


la matrice de découplage correspondante sont obtenus comme : 


0.7278(20.463s + 1) 
(17.938s+1) 


D(s) = 0.352e 01595 1 —0.0028 
23.238(18.29s+1) _,,,0, —26.457(4.128s +1) 
mt F4 e 


(7.27195+1) ° (14.521 + 1) 


1 —1437s 0. 0048e 22035 


—6.4305s 1 


Les éléments de la matrice RGA et la valeur de l'indice NI sont : 


1.96 —0.66 —0.30 
RGA =|-—0.67 1.89 —0.22| et NI = 0.4006 
—0.29 —0.23 1.52 


L'appariement recommandé est y. avec u:,y, avec u, et y; avec u;,. parmi les six configurations 
possibles. 

Dans cet exemple, nous sommes appelés à faire une comparaison des méthodes BLT (Biggest Log- 
modulusTuning) [106],Lee et al. (2004)[107] et CMAES-FOPI (2014) [108], qui sont les méthodes les 


plus populaire pour le réglage du contrôleur PI décentralisé. 
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Les paramètres optimisés des contrôleurs fractionnaires FO-PI obtenus par la méthode proposée sont 
organisés respectivement dans le Tableau 3.6. 


Tableau 3-6:Paramètres optimaux du régulateur FO-PI décentralisé ainsi que ceux du 
méthodes réglage précédentes. 


MANSRES Paramètres du contrôleur 
de réglages 
critère de 
OR Column Kp1 K p2 KP3 Ki K 2 K;3 Li L2 L3 performance 
IAE 

Proposed 

PSO-FOPI 0.00017 -0.093 4.99 0.195 -0.112 2.407 0.9997 0.9758 0.998 48.3268 
controller 

BLT 

(Luyben, 1.510 -0.295 2.63 0.092 -0.016 0.397 - - - 343.76 
1986)[106] 

Lee et al. 
(2004)[107] 0.593 -0.124 3.22 0.172 -0.043 0.42 - - - 280.36 
CMAES-FOPI 
(2014)[108] -0.3663 -0.6279 4.999 0.2178 -0.1456 4.9998 0.92061 0.9886 0.7971 61.0379 


La figure 3.12 (a,b,c)montre la réponse indicielle des première, deuxième et troisième sorties 


(V1, Y2etY;) du système contrôlé avec des contrôleurs fractionnaires et classiques décentralisés. 


La modélisation et la simulation de processus multivariables de grandes dimensions sont très 
difficiles ou peu pratiques en raison des fortes interactions entre les entrés et Les sorties du processus. 
Les chercheurs ont toujours eu recours à l'élimination des effets de ces interactions en concevant des 
dispositifs de découplage adaptés, transformant ainsi le processus en boucles indépendantes. Dans 


notre recherche présentée, cela a été réalisé, et nous pouvons résumer nos résultats comme suit : 


% La technique de réglage des paramètres des contrôleurs PSO-FOPI basée sur la méthode 
d’optimisatio PSO est plus supérieure que les méthodes mentionnées ci-dessus sont obtenues à 


partir de la littérature existante. 


“ Des changements de consigne ont été introduits séquentiellement dans les boucles 
individuelles (Une variation de consigne de 0 à 1, est appliquée sur rl, à l’instant t-0s. Une 
variation de consigne de 0 à 1, est appliquée sur r2, à l’instant t-200 s. Une variation de 


consigne de 0 à 1, est appliquée sur r3, à l’instant t-400 s. ) 


“ La performance supérieure du contrôleur proposé est apparente dans la Figure 3-11et le 


Tableau 6. 


“ Les résultats des réponses de sortie du controller PSO-FOPI proposé avec le critère IAE était 
plus favorable par rapport aux régulateurs PI conventionnels. 

L'efficacité et l'applicabilité des contrôleurs SPO-FOPI proposés pour la colonne de distillation 

Ogunnaike et Ray (OR) (procédé multivaariable complexe de trois entrées et trois sorties) ont été 


démontrées, bien que la mise en œuvre soit difficile. 
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Figure 3-12:Les Réponses indicielle du système contrôlé: (a) une consigne de 1 est appliquée 
sur Ÿ,, a l’instant t-0 s. (b) une consigne de 1 est appliquée sur Ÿ,, a l’instant t-100s. (c) une 


consigne de 1 est appliquée sur Y3, a l’instant t-200 s. 
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3.5 Conclusion 


Dans ce chapitre, nous avons résolu le problème de réglage optimal des paramètres de 
contrôleur PI d'ordre fractionnaire pour les systèmes MIMO en utilisant l'algorithme PSO. Les 
interactions ont été minimisées par la procédure de conception proposée du contrôleur avec de 
découplage. La conception du contrôleur PSO-FOPID a été étudiée en l'appliquant aux processus de 
colonne de distillation (Type: processus 2x2 et 3x3 ) pour obtenir les performances souhaitées . Les 
résultats de simulation de la méthode proposée et de diverses méthodes (NDT-PI and SIMC-PI } sont 
comparés par la vérification des normes de performance du contrôleur. Cette méthode donne un 
meilleur suivi du point de consigne et rejette rapidement les perturbations par des interactions 
minimales dans les boucles individuelles du système MIMO. On peut dire que les performances du 
régulateur PI d'ordre fractionnaire basé sur l'approche diffusive sont nettement supérieures aux 


performances des régulateurs PI classiques. 


Chapitre 4: Synthèse des contrôleurs fractionnaires pour les 
systèmes multivariables basée sur la méthode de commande 
IMC. 


4.1 Introduction: 


Généralement la synthèse de lois de commande pour les systèmes d'ordre élevé mène à 
des contrôleurs plus complexes, ce qui n'est pas souhaité en domaine pratique. Car concevoir 
de tels contrôleurs est en réalité très coûteux du point de vue économique. Alors, pour 
résoudre ce problème, on a recours à des méthodes de simplification des systèmes d'ordre 
élevé en structures simples (modèles simples) et avec un ordre réduit, à travers lesquelles il 
est possible synthétiser les contrôleurs. Et dans ce contexte, le choix d'une méthode 
appropriée pour synthétiser les contrôleurs est très important pour obtenir les performances 
désirées et assurer la stabilité et la robustesse de la commande. Le choix d'une structure de 
commande par modèle interne (IMC) est l'une des meilleures structures de commande qui 
permettent d'atteindre ces objectifs. 

Récemment, l'attention est accordée à la synthèse d'une loi de commande d'ordre non 
entier basé sur la structure de commande par modèle interne (IMC) pour les systèmes 
monovariables, par exemple Bettayeb et Mansouri [109-11] obtient un nouveau type de 
contrôleurs appelé " contrôleur PID classique avec un filtre d'ordre fractionnaire " (en anglais 
appelé "FOF-PID Controller") [12-110]. 

Dans ce chapitre, nous présentons le principe de la commande par modèle interne ainsi que 
la méthode de synthèse du contrôleur pour la commande des systèmes monovariable stable. 
Nous proposons une structure de commande par modèle interne (IMC) qui permet la synthèse 
de contrôleur d'ordre fractionnaire pour les systèmes multivariable. La méthode proposée est 


appliquée au système TITO avec decoupleur inversé. 


4.2 Synthèse par l'approche de la commande à modèle interne (IMC) 
pour les systèmes SISO 


4.2.1 Commande par modèle interne: 


Ces dernières années, la méthode de commande par modèle interne (IMC) a acquise un 
intérêt croissant, en raison de la simplicité de synthèse des contrôleurs et la robustesse 
inhérente à sa structure qu'elle offre [112-113]. La simplicité de cette approche réside dans le 


fait qu'il n'existe qu'un seul paramètre à régler, directement, lie à la constante de temps en 
q q P 8 p 
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boucle fermée définie par l'utilisateur. La méthode par IMC regroupe les avantage de la 
boucle ouverte qui se présentent par la facilite de la synthèse du contrôleur et la possibilité de 
tenir compte de la robustesse et de pouvoir traiter les procèdes avec retard, et les avantages de 
la boucle fermée qui sont la possibilité d'obtenir une erreur statique nulle en régime 
permanent avec des consignes de type échelon en présence de perturbations ainsi que des 


erreurs de modélisations du procède à contrôleur 


Contrôleur Système 


Modèle 


Figure 4-1:Schéma de commande par modèle interne (IMC) 


La figure 4-1 représente la structure de commande par modèle interne, où G(s) représente 
la fonction de transfert du système à commander , et son modèle est donné par la fonction de 
transfert Gn(s). Dans cette structure de commande r(s) représente le signal de référence 
qui est comparé à la l'écart £(s) entre la sortie réelle y(s) et la sortie mesurée où ym(s): 
E(s) = y(s) — ym(s), la différence est introduite à l'entrée du contrôleur C;,,.(s) dont la 
sortie est l'entrée du système g(s) et de son modèle gm(s) ce qui permet à ce schéma de 
prendre en compte des erreurs de modélisation, la grandeur de perturbation d(s) est ajoutée à 
la sortie du procédé. 

Remarque : 
Ainsi, dans le cas parfait où (d(s) = 0) c.-à-d. le système n'est pas affecté par les 
perturbations. Et aussi e(s) = y(s) — ym(s) = 0 c.-à-d. que la sortie du modèle G,;,(s) 
représente parfaitement le système G(s)cette structure de commande devient simplement une 
structure de commande en boucle ouverte (le contrôleur fonctionne en boucle ouverte). 
D'après le schéma de la Figure 4-2, la sortie y(s) peut être exprimée en fonction de la 
perturbations d(s) et la référence r(s) par : 

G(S)Cime(S) 1 — Cime(S)Gm(s) 
ETPOILOE "OUEST TO ET) 


La Figure 4-2 représente la commande classique à retour unitaire. 


y(s) = d(s) (4.1) 
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Système 


&(s) 


Figure 4-2:Schéma de commande classique à retour unitaire. 


La sortie y(s) dans le cas d'une commande en boucle fermée classique illustrée à la Figure 4- 
2 est donnée par 


G(s)C(s) 


+600)" + 


y(s) = d(s) (4.2) 


1 + G(s)C(s) 
On peut définir la relation entre le contrôleur classique C(s) et le contrôleur à modèle interne 


Cimc(S) à travers les relations (4.1) et (4.2). 


Cime(S) 
61715, 60m) ” — 
ee (2) 


1 +Gn(s)C(s) 
Les relations de l'équation (4.3) permettent de passer de la commande à modèle interne à la 
commande classique et vise et versa. Par conséquent, les performances qui peuvent être 
obtenues en utilisant un contrôleur à modèle interne ne sont, a priori, ni moins bonne, ni 
meilleures que celles obtenues en utilisant un contrôleur classique, puisque ces deux structure 


sont équivalentes. 


4.2.2 Synthèse du contrôleur IMC pour le système stable: 


Si nous supposons que G,(s) = G(s) signifie que le modèle représente parfaitement le 


système, alors l'équation (4.1) devient 


Ys) = G(S)Cimets) T(S) + (1 — Cimc(S)Gm(S)) d(s) (4.4) 
La synthèse du contrôleur IMC se réalise en deux étapes [112-113]: 
Etape 1 : La synthèse du contrôleur C;,.(s) est simple, elle utilise le principe de la 
commande en boucle ouverte. Néanmoins, à cause de l'inversion du modèle, lorsque G; (s) 
possède un retard ou des zéros instables il est factorisé en deux parties inversible et non 


inversible selon : 


Gm(s) = Gm (S)Gm (s) (4.5) 
Où: GX (s) est la partie non inversible de GA (s), et inclut les zéros instables et l'élément 
retard. En plus le gain statique de ce transfert doit être égal à 1 (G} (0) = 1). Le reste de la 
fonction de transfert G,(s), est G; (s) qui s'appelle la partie inversible . C'est cette partie qui 
est utilisée pour calculer le contrôleur Cyn.(S). 


Etape 2 : la fonction de transfert du contrôleur Cimc(s) est donné par 
1 
Gm (S) 


Le filtre f(s) représente la fonction de transfert désirée du système en boucle fermée, appelé 


Cime(S) = fs) (4.6) 


aussi "filtre de robustesse" par Morari [112-113]. Il est généralement choisi égal à : 


1 
f(s) = G+ry (4.7) 


Où: +. étant la constante de temps qui définit la durée du régime transitoire, et le paramètre r est un 
nombre entier positif, il est choisi de sorte que le contrôleur finalement obtenu soit propre et 
causal. 

En utilisant l'équation (4.6), le contrôleur à modèle interne C;m.(s) est calculé et en substituant 


dans l'équation (4.1) le contrôleur classique C(s) peut être retrouvé. 


4.2.3 Synthèse des contrôleurs IMC-PID-FOF/ IMC-FOPID-FOF 


Afin d'imposer un comportement fractionnaire au système en boucle fermée, comme 
suggéré par Bettayeb et Mansouri dans leurs travaux récents [114-115-11-12]. En choisissant 
la fonction idéale de Bode en boucle fermée donnée par l'équation (1.29) comme filtre de 
robustesse lors de l'application de la méthode IMC. 

Dans ce cas, et afin d'avoir en boucle fermée un comportement fractionnaire robuste, On 
choisit la fonction de transfert de filtre f(s) d'ordre fractionnaire égal à la fonction idéale de 
Bode en boucle fermée. Le contrôleur résultant est appelé un IMC-PID-FOF , qui est 
l'abréviation du correcteur PID classique en série avec un filtre d'ordre fractionnaire (FOF) 

La méthode de synthèse du contrôleur IMC-PID-FOF proposée se résume en trois étapes 
principales [11]: 

Etape 1 : 

On choisit la fonction de transfert de filtre f(s) égal à la fonction idéale de Bode en boucle 


fermée comme modèle de référence lors de la synthèse du contrôleur : 


fs) = Grer(s) = 0<a<1 (4.8) 


1+7.41 


Ainsi si le modèle G,,(s) représente parfaitement le système G(s) Il est factoriser en deux 
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parties inversées et non réversibles, selon l'équation (4.5) avec G (0) = 1 pour s'assurer qu'il 
n'y a pas d'erreur constante. 

Etape 2 : 

En utilisant l'équation (4.6), la fonction de transfert du contrôleur Cync(s) est calculé et donné 
par : 


1 1 


1 
Cimc(s) = OUS aOirrnt (4.9) 


Les deux paramètres T. et « sont choisis afin de spécifier respectivement une marge de phase 
Pm (qui est liée au dépassement) et la fréquence au gain unité w. (qui est liée au temps de 


réponse) du f(s), respectivement [11-116]: 


= Te 1 and Te = (4.10) 
Etape 3 : 
La fonction de transfert du contrôleur équivalent C(s) est donnée par 
C(s) = | (4.11) 


Gn (SJ + Tet1 — Gr (s)] 
e Dans le cas du modèle G,,(s) à ordre entier, l'expression analytique du contrôleur C(s) peut 
être mise sous forme d'un contrôleur PID d'ordre entier en cascade avec un filtre d'ordre 


fractionnaire (FOF), d'où la nomination du contrôleur IMC-PID-FOF [11]. 


re 


1 
(s) kp Gi + Ts + Tas) (4.12) 
FPT 


Conventional PID 


e Dans le cas du modèle G,,(s) à ordre non entier, l'expression analytique du contrôleur 
C(s) peut être mise sous forme d'un contrôleur PID d'ordre fractionnaire (FOPID) en 
cascade avec un filtre d'ordre fractionnaire (FOF), d'où la nomination du contrôleur 


IMC-FOPID-FOF [11]. 


1 
C(s) = m(s) kp (: ar Tast) (4.13) 
FOF Ts 


PID Fractionnaire 
4.2.3.1 Cas d'un système décrit par un modèle rationnel du premier ordre à retard (FOPDT): 


Considérons maintenant le modèle rationnel (d'ordre entier) du premier ordre à retard (En 
anglais: The First Order Plus Dead Time (FOPDT) model) donné par la fonction du transfert 


suivant: 


K 
Gm (S) = : ee pi (4.14) 
m 


Le contrôleur IMC est donne par: 


1ETS 1 


Cime(s) = Kn, Tel (4.15) 
Le contrôleur équivalent C(s) est donné par: 
1+Tns 
— (4.16) 


C = —————————————————— 
(s) Kn[1 + Ter = e-tms] 


Pour approximer le retard, deux méthodes d'approximation peuvent être utilisées. La première 
est l'approximation de Taylor de premier ordre qui permet d'avoir un contrôleur IMC-FOF-PI 
et la deuxième est l'approximation de Padé de premier ordre qui aboutit à un contrôleur IMC- 
FOF-PID 
e Approximation du retard par la formule de Taylor du premier ordre 
e7êms = 1 — Oms (4.17) 
En substituant le retard approximé dans l'équation (4.16), la fonction de transfert du 


contrôleur équivalent C(s) est donnée par: 
1+TnS 


Ce — 
6 Km + T.2+1 + Oms] 


(4.18) 


On note qu'il peut être mise sous forme d'un contrôleur PID d'ordre entier en série avec un 


filtre d'ordre fractionnaire (FOF) (appelé : contrôleur IMC-PI-FOF) : 


C(s) 1 1 Ç + s) (419) 
S)=— —— : 
1+(T/0m}s® K,0m s 
filtre fractionnaire contrôleur PI classique 
FOF PI 
Qui peut s'écrire comme: 
C(s) : ke(1+=) (4.20) 
S = ———— C — : 
1 +75 T;s 
———— = à 
filtre fractionnaire contrôleur PI classique 
FOF PI 
où: 
Tn 
K. _— 
5 K,Om 
T=T, (4.21) 
Tee T./0m 


e Approximation du retard par la formule de Padé du premier ordre 


Si on approximé le retard au selon l'approximation de Padé du premier ordre : 


om 
-êms I ZE 4.22 
e = 2 (4.22) 

1+5-sS 


Le contrôleur Cy.(s) devient : 


Om 
A+7ns) (+) 
Km (A+t.t+1) 


Cime(S) = (4:23) 


Et le contrôleur C(s) correspondant est : 


(1+TnsS) (1 4 Ms) 
s) = don (4.24) 
Kn [re + Ts] 
Encore une fois, On note qu'il peut être mise sous forme d'un contrôleur PID d'ordre entier en série 


avec un filtre d'ordre fractionnaire (FOF) (appelé : contrôleur IMC-PID-FOF) : 


Éeÿ = 1 2Tn + 0m : 1 " Tn0m re 
AA NT Pis Kn0m 2Tn + Om 2T,, +0m° os 
1 + 0 S 
m 2 
filtre fractionnaire contrôleur PID classique 
FOF PID 
qui peut s'écrire comme 
C(sS) = 1 Ke(1++T ) 4.26 
: _ 1+ TFs® . T;s as ( ù ) 
Rene Pré EEE DRE 
filtre fractionnaire contrôleur PID classique 
FOF PID 
où: 
__ 2Tn + 0m 
SO Knom 
2Tn + 0m 
T, = LR 
2 
(4.27) 
: Tn0m 
47 27, +0m 
2Te 
= Om 


4.2.3.2 Cas d'un système décrit par un modèle d’ordre non entier à retard de type I (NIOPTD- 
D : 


Considérons un modèle d'ordre non entier de premier ordre plus un temps de retard 


(NIOPTD-D), décrit par la fonction de transfert : 


Km 


x = —On 
Gu(s) NIOPTD-I — Tas+1° À 0<y <2 (4.28) 


Le contrôleur IMC calculé par l'équation (4.6) est donné par : 
1+TnS* 1 


Cimc(S) = Rs Leri (4.29) 
Le contrôleur équivalent C(s) est donné par: 
1+TyhsS’ 
C(s) = = — (4.30) 


m li + Teat1 _ e-68ms] 


e Approximation du retard par la formule de Taylor du premier ordre 
e7Êms = 1 — Oms (4.31) 
En substituant le retard approximé dans l'équation (4.17), la fonction de transfert du 
contrôleur équivalent C(s) est donnée par: 


C S) — 


—————— 4,32 
mlec2+1 + 0ms] ( ) 


On note qu'il peut être mise sous forme d'un contrôleur PI d'ordre fractionnaire (FOPT) en 


série avec un filtre d'ordre fractionnaire (FOF), (appelé : contrôleur IMC-FOPI-FOF) : 


C(s) = 


sr” 1 1+TmS) 
(7) (4.33) 


1+(T/0m}s®t K,0m sY 


Len” | LE” | 
filtre fractionnaire contrôleur PI classique 
FOF PI 


Qui s'écrit finalement sous la forme générale: 


sY-1 


C(s) Kc Gi + (4.34) 


A + Ts 


a” a ———— 
filtre fractionnaire PI d'ordre fractionnaire 
FOF FOPI 


(4.35) 


Tr= T./0m 


| Système d'ordre élevé 


Première Démarche L L Deuxième Démarche 


Réduction du modèle original | Réduction du modéle original par 


par un modèle d'ordre entier un modéle d'ordre fractionnaire 
de la forme FOPDT de type I (NIOPDTI-I) 


Calcul du contrôleur d'ordre Calcul du contrôleur d'ordre 
fractionnaire fractionnaire 
Mise sous la forme Mise sous la forme 
IMC-FOF-PID IMC-FOF-FOPID 


Figure 4-3: Démarches à suivre pour la synthèse des contrôleurs fractionnaires pour les 


systèmes d'ordre élevé 


4.2.4 Exemple d'illustration. 
Le système d'ordre entier supérieur suivant est utilisé dans cette simulation [117]: 


1 


GS) = 
5) (s + 1)(0.2s + 1)(0.04s + 1)(0.008s + 1)) 
Le modèle NIOPDT-I pour ce système est obtenu comme: 
0.99932 
* — —0.1922s 
Cu) = Losazsor +1 


Les réponses indicielles du système d'ordre entier supérieur G(s) et le modèle NIOPDT-I 


sont illustrés à la Figure 4-4. 


z 
(s + 1)(0.2s + 1)(0.04s + 1)(0.008s + 1)) 


——G(Cs) 


Sortie 


0.099932 —02922= 


—— CG) 6822510157 +1 


1 1 1 1 1 1 n 1 1 
e} 2 4 6 8 10 12 14 16 18 20 
Time (s) 


Figure 4-4:Réponses indicielles du système d'ordre entier supérieur et le modèle NIOPDT-I 


Les spécifications de conception du domaine fréquentiel du système en boucle fermée sont 
choisies comme nn, = 80 et wc = 0.4rad/s. 

En utilisant l'équation. (4.10), les paramètres de la fonction de transfert idéale du Bode sont 
trouvés comme: &æ = 0.1111et 7, = 2.7679 et donc l'équation. (4.33) qui exprime le 


contrôleur proposé IMC FOPI-FOF a la forme suivante : 


C{s) = 50182 1.0842 (1 < 1 ) 
LecT” (2.7679/0.1922)50:1111 0.99932 + 0.1922 1.084251-0212 
filtre fractionnaire PI d'ordre fractionnaire 
FOF FOPI 


Les figures 4-5 et 4-6 montrent la réponse indicielle du système commandé par IMC FOPI-FOF 


Controller pour différentes valeurs @,, et wc des paramètres du modèle de référence. 


Réponse Indicielle 


We= 0.1 rad/s 
We= 0.2 rad/s 


We= 0.3 rad/s 


Sortie 


We= 0.4 rad/s 


We= 0.5 rad/s 


O 20 40 60 80 100 
Temps (sec) 


Figure 4-5:Réponse indicielle du système commandé pour différentes valeurs du paramètre 


We. 


Réponse Indicielle 


4 
————— W,.= 0.4 rad/s et 
© 0.8 : ?m 
= —— W,= 0.4 rad/s et p,, 
a 0.6 = "# = 0.4 rad/s et Pm 7 
0.4 —Ww,= 0.4 rad/s et Pm . ] 
0.2 L 


O 20 40 60 80 100 
Temps (sec) 


Figure 4-6:Réponse indicielle du système commandé pour différentes valeurs du paramètre 


Pm- 


4.3 Synthèse par l'approche de la commande IMC pour les systèmes 
(MIMO) 


Dans cette partie, la structure de commande basée sur l’approche IMC pour les systèmes 
SISO sera généralisée à la commande des systèmes MIMO afin d'améliorer les performances 
du système en boucle fermée. L'effet de toutes les boucles est considéré comme une 
perturbation sur la boucle de commande considérée, d'où l'idée d'utiliser découplage , pour 
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séparer les interactions entre les variables du système MIMO. 


4.3.1 La commande IMC avec découplage inversé 


Dans notre méthode proposée, nous utilisons la technique de découplage inverse, où le 
système générale (système et le découplage) est directement transformé et devient un 
système diagonal. Dans ce cas, les techniques de synthèse des contrôleurs monovariables 
peuvent être adaptées pour la commande des sous-systèmes obtenus indépendamment l'un 
de l'autre. 

La contribution de cette étude est d'imposer un comportement fractionnaire (commande 
d'ordre fractionnaire) pour les boucles de commande du système en boucle fermée. Cela se 
fait, en imposant le modèle de référence, pour chaque boucle, comme étant la fonction idéale 
de Bode en boucle fermée pour synthétiser le contrôleur IMC (IMC: Internal Model 
controler). Par la suite, le contrôleur équivalent dans la structure de commande multi-boucle 
standard, est déterminé. 

Maintenant, la méthode de synthèse pour le contrôleur IMC-FOPI-FOF[16], présentée 
dans la section (4.2) pour la commande des systèmes monovariable est généralisée à chaque 
boucle des boucles séparées. En fin, on obtient alors un contrôleur diagonal pour le système 
multivariable (dite: contrôleur multi-boucle d'ordre fractionnaire ) 

La structure de commande multi-boucle pour un système multivariable ( MIMO: de deux 


entrées et de deux sorties 2x2 ) par modèle interne est illustrée à la Figure 4-7. 


Où: Le système à contrôler et son découplage inversé sont représentés respectivement par 
G(s) et D(s). Les systèmes découplés sont représentés par Q(s), et leurs modèles équivalents 
sont Gu(s). 

La méthode de la synthèse du contrôleur multi-boucle fractionnaire IMC-FOPID-FOF proposé pour 
un système MIMO (de deux entrées et de deux sorties, en anglais: TITO system) se réalise en 


trois étapes : 
4.3.1.1 Etape 1: Modélisation des systèmes découplée par un modèle de type I (NIOPTD-T): 


Nous concevons une technique de découplage inversé pour décomposer un processus 
multivariable stable en une série de sous-systèmes mono-boucle indépendants. La conception 
d’une matrice de transfert (s), tel que D(s) est en série avec G(s) produisant une matrice de 
transfert Q(s) diagonale. 


q11(s) 0 | 
0 q22(s) 


Les systèmes découplés Q(s) = diag(qa1{s),q22(s)) sont des systèmes d'ordre élevé qui sont 


Q(s) = G(s).D(s) = | (4.36) 


réduits à des modèles d'ordre fractionnaire de type I (en anglais: NIPODT-D , a été exprimée 
dans l'équation (4.28) ci-dessus, où: 


Gn (s) 11 | 0 


d11(S) 0 
GS) = G(s).D =| ]= . 4.37 
m(s) “DE 0 q22(s) | 0 Gn(S)22 ( ) 
avec : 
. Km ii un. 
Gu(S) vropro-1 = Tausü+1 on LS (4.38) 


Si q11(S)= G11(5)* et gq22(5)= G22(s)* signifie que les modèles représentent parfaitement les 
systèmes découplés. 

Dans l'étape de modélisation, les paramètres de chaque modèle sont donnés par le vecteur: 

X1 = [Kon 1 Im 1 0mMm: Yu pour le premier modèle Gy(s)';;, (4.39) 

X2 = [Kn22 TM 22 O0Ma2 Y,,] Pour le deuxième modèle Gy(s) ;; 

Les conditions de ces paramètres optimisés sont données par : 


Kiun < Km 1,2 < Ke 
0 <Tm 1,2 < Las 
0 < Om: 2 < On max 
0<Yy,, <2 


(4.40) 


La fonction objectif pour la réduction de modèle basée sur l'utilisation de l'algorithme 


d'optimisation PSO peut être définie comme: 


n 
1 
Le (Kn 11 TM 31, OM, Vu ) = => On — ÿ,)? pour le premier modèle . (4.41) 
i=1 


n 
1 
Jmin (Kn 22, Tm 22, 0M22,Y 2 ) = = ù O2 — ÿ,)? Pour le deuxième modèle. (4.42) 
i=1 


Avec: 

Y1 = us); Y2 = q22(s); TION Ÿ2 = Gu(S)' >» 
La méthode de la modélisation des systèmes découplés avec des modèles d’ordre fractionnaire 
équivalents de type I (NIOPDT-T) peut être représentée sur la Figure 4-8. (Pour plus de détails sur la 


modélisation du systèmes, voir [16]). 


Z - 
RE > Va 
| Joraies CH 22: T978 12, Onrtss, Vas D = > Gr 2 H——- 
= ++. 
= S E 4 TITO systerm 
— Découplasge inversé 
| æ BP | 


nr 
4 


2? D 
pes) MR GS) — diag(qs3(5),422 C5) 
: 7 Jin (Km22,: Tom 22: O2, V3 > = . > G2 — 227 | = - 
ti 
= 


n— 
— 


Figure 4-8:Modélisation des systèmes découplée à modèles d’ordre non entier à retard de 


type I (NIOPTD-I) . 


Le procédé multivariable peut être contrôlé à l'aide des contrôleurs de boucles indépendants 


appelés: contrôleurs décentralisés (ou: Le contrôleurs multi-boucle) 


4.3.1.2 Etape 2: Synthèse du contrôleur multi-boucle fractionnaire IMC-FOPI-FOF 


Le contrôleur multi-boucle IMC est donné par la matrice de transfert: 


Cime(S) = diaglCime(S)11, Cimc(S)22] : Es C; to. 


Afin de calculer les contrôleurs IMC monovariables Cimc(s)i;, i = 1,2 on procède d'abord à 


(4.43) 


factoriser chaque élément G,, (s);; en deux parties : inversible et non inversible [49-118-1191. 
Gm (S)ä = Em (S)i Gm'(S)à avec: Gm'(0)ÿ = 1 (4.44) 


Chaque contrôleur de la matrice de transfert (4.43) est alors calculé comme suit : 
1 
Cime(Si = fi | i= 1,2 (4.45) 
me Ju Gn (s)i: r [22 


Avec: f(s);; est le modèle de référence imposé à chaque boucle (u;, y;). Dans notre étude, on 
propose d'imposer un comportement fractionnaire pour chaque boucle en choisissant la 


fonction idéale de Bode en boucle fermée comme modèle de référence [118-49]. 


1 
1 + Tatüti / 


fi = 0O<a;<1 (4.46) 


Comme nous l'avons fait précédemment dans le cas des systèmes monovariables. Les deux 
paramètres Tai et æi sont choisis alors en spéciant la marge de phase et la fréquence de 
coupure désirées (qui est liée au dépassement et temps de réponse respectivement) pour 


chaque boucle en boucle fermée : 


T — ini 1 
Œii — r/2. — 1 and Toii — PPNCTEE (4. 47) 
Le contrôleur IMC calculé par l'équation (4.45) est donné par : 
1 + GrniiS'À 1 
Cimc(S)i = T— (4.48) 


Ki 1+Toutiiti 
4.3.1.3 Etape 3: l'implémentation 


La figure 4-9 illustré la structure général de commande multi-boucle classique où G(s) est 
le procédé multivariable à contrôler, il est donné par l'équation (2.1), D(s) est le découplage 
inversé il est donné par l'équation (2.1), et C(s) est le contrôleur multi-boucle d'ordre 


fractionnaire à déterminer, il est donné par l'équation (4.49) : 


T2 


M Découplage inversé 
D(s) 


Figure 4-9:Structure général de la commande multi-boucle classique. 


(4.49) 


C(s) = diag[C(s)11,C(S)22] = Ee | 


L'équation (4.49) obtenue est le résultat de l'équivalence entre les deux structures de 


commande IMC multi -boucle et classique respectivement des figures (4.7) et (4.9), 


Cimc(S )i 


CS = —— (4.50) 
. 1 Cimc(S)üiGm(S); 
Le contrôleur équivalent C(s);; est donné par : 
1 + Gi is" Ù 
C(S}i = — (4.51) 


Kou LL + Tete 6mits] 


Pour le retard e-%"iS, on utilise une approximation de Taylor au premier ordre, et la fonction de 
pp y P 


transfert du contrôleur C(s);; devient la suivante 


ii—1 ii 
Re (4.52) 
M 1+(rx/0mii)s KrOmii sr | 
filtre fractionnaire contrôleur PI classique 
FOF PI 


Le contrôleur multi-boucle à déterminer mise sous forme d'un contrôleur PI d'ordre 
fractionnaire (FOPI) en série avec un filtre d'ordre fractionnaire (FOF), (appelé : contrôleur 


multi-boucle fractionnaire IMC-FOPI-FOF). 


4.4 Simulation et Application 


Un exemple est présenté afin d'illustrer la méthode de synthèse du contrôleur multi-boucle 
proposé pour la commande du système multivariable (système de deux en entrées et deux 
sorties). Les performances des contrôleurs obtenus sont comparées à celles des autres 


contrôleurs existants dans la littérature. 


128 _, 189 

= 6€ € 

_ |16.75+1 215+1 
AIS 6e … 14 


1095+1°  1445+1° 
La matrice de gain relatif RGA et la valeur de l'indice Niederlinski Ni sont obtenues à l'aide 
des équations (2.10 et 2.20) respectivement : 


2.0094 —1.0094 


RGA(G(0)) = —1.0094 2.0094 


| NI-0.4977 


La configuration de commande appropriée est : u: contrôle y. et u contrôle y2. 
En appliquant la technique du découplage inversé pour éliminer les interactions du système. 


La matrice correspondante est obtenue comme suit : 
1.477(16.75+1) .. 
en ne à 


1 
Di) = 03444541) nn 
7 109511 : 
Les systèmes découplées de Q(s) = G(s)2x2 X D(s)2x2 = diag(q11(5), q22(S)) sont : 
1 12.8 6.43(14.45+1) 
PE EC CES IE EEE VS CE VEUT) 


(10.95 + 1)(215 +1) 


1 19.4 9.7453(16.75s + 1) w) 


d22(s) = 2 e 
0.5017(14.4s + 1)(16.75+1)|\ (14.45+1) (215 + 1)(10.9s +1) 
(10.9s + 1)(21s + 1) 


Après l'étape de modélisation proposée du système découplé, les modèles d'ordre 


fractionnaire de type NIOPDT-I ont été obtenus comme suit : 


13.145 
x = —1.2107 
Emi SE LE osggsos80 1 
_19.7716 
* = —3.1034 
Gmez SE Lies 1 


Les réponses indicielles sont données sur la Figure.4.10. Pour évaluer l'affinité entre le modèle d'ordre 


réduit (modèle d'ordre fractionnaire proposé) et chaque sous-système découplé équivalent. 


Réponse indicielle 


_—-—- Original @, 1) 


6 Fractional order model (NIOPTD-I nn) + 
4 = 
2 = 
O 1 i 1 i 
O 20 A0 60 80 100 
Temps (sec) 
(a) 
ü Réponse indicielle 


_--—= Original (ad) 
Fractional order model (NIOPTD-1,5) 


O 20 AO 60 80 100 
Temps (sec) 


(b) 
Figure 4-10:Réponse indicielle en boucle ouverte de la fonctions de transfert du chaque système 


découplée réel et de son model approximation (NIOPTD-T). 


Pour cet exemple de simulation, si les valeurs souhaitées du modèle de référence en termes 

de fréquences de coupure et de marge de phase sont respectivement: 

Pm11 = Pm22 = 81° Et Woy1 = Wo22 = 0.4 rad/s. 
Et en utilisant les équations (4.46 et 4.47), nous avons trouvé les valeurs des constantes de 
temps : Te11—Te22— 2.7399 et les valeurs de l'ordre fractionnaire &11-æ22 sont 0,10. Ces valeurs 
sont choisies pour répondre aux spécifications souhaitées du modèle de référence en boucle 
fermée dans le domaine fréquentiel. Le contrôleur IMC FOPI-FOF obtenu avec cette 
méthode proposée est présenté dans le tableau 4.1. 

Dans ce travail, la méthode proposée est comparée aux derniers travaux avancés et aux 
nouveaux résultats obtenus par les chercheurs dans le domaine du contrôle des systèmes 
multivariables tels que les techniques proposées par Chuong et al.[120][75] , et la méthode 
proposée par Chananchai et al.[121]. Les contrôleurs conçus par ces techniques dans [120, 


75, 121] pour le même exemple de simulation sont répertoriés dans le tableau 4.1. 


Tableau 4-1:Les contrôleurs multi-boucles résultants avec différentes méthodes de réglage. 


Méthode de réglage Les contrôleurs résultants 
Le contrôleur multi-boucle HE 1 
; CO = T7. 263s010 1 0027 (1 ET] PE) 
IMC proposé 
57001103 _— ” 1 
DE Le d CCSd © oran (C0-28) (1 + 5 gs) 


Fractional-SDSP 


1 1 
C(S)11 = 0.2882 5e (1 a 
Chuong et al. 2019 [120]. 


0.1286 7.1079508714 


1 1 
C(S)22 = —0.1212 ess (1 Faso + 0.6385509%8) PS 
ical- 1 1 
Classical-SDSP C(s)22 = 2.0205 (1 + 19.5275s +2. 16175) 9.192s +1 
Chuong et al. 2019 [75]. ‘ l 
C(S)22 = —0.330 (1+ : + 0.14285) : 
ARE 14.481157 © ‘/8.2745+1 
Classical-PID C(s) (0.460952 + 0.5387s + 0.0774) 
S) 44 = 
Chananchai et al. 2018 [121]. à 
(—0.12135? — 0.1583s — 0.0384) 
s 
: 1 1 
with forward controllers : G,(S)11 = an ; Gp(S)22 = en 


Les résultats de la réponse indicielle du système en boucle fermée des deux boucles (bouclel et 


boucle2) est illustré sur la Figure 4-11 (a et b), respectivement, avec un changement du point de 
consigne, les changements du point de consigne sont séquentiels, pour la boucle 1 à t-0 sec et la 


boucle 2 à t-100 sec. 
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Figure 4-11:Réponse indicielle du système TITO en boucle fermée avec différentes méthodes 


de réglage. 


La Figure 4-12 (a et b), montre les signaux du contrôleur multi-boucle proposé appliqués 


au système commandé. 


signal de commande, 


0.8 ï. T T T T T T T T T 
k 
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Figure 4-12:Les signaux de commande appliquées. 


Afin d'évaluer les performances de la réponse de système en boucle fermée dans le 
domaine temporel (temps de stabilisation, dépassements, L'intégrale de la valeur absolue de 


l'erreur (IAE)) par les quatre méthodes, nous avons résumé dans le Tableau 4.2. 


Tableau 4-2:Résultats de l'analyse des performances du système en boucle fermée dans le 


domaine temporel. 


Méthode de réglage Temps de réponse (sec) Dépassement ( %) IAE 
Boucle; Boucle; Boucle; Boucle; Boucle; + Boucle; 
1-71) | OZ-r2) | (y1-r1) | (2-r2) 
Contrôleur IMC proposé 9.98 16.77 7.6 10.01 9.053 
Fractional-SDSP 14.86 19.99 3.11 5.07 9.221 
(Chuong et al. 2019) 
Classical-SDSP 26.14 25.32 13.09 15.89 10.963 
(Chuong et al. 2019) 
Classical-PID 14.92 40.82 0.52 41.03 12.08 
(Chananchai et al. 2018) 


Pour démontrer l'efficacité de la méthode de conception proposée, les réponses indicielles 


de la simulation sont présentées pour faire une comparaison des performances du système en 


boucle fermée dans le domaine temporel avec les quatre méthodes de réglage. 


Nous pouvons voir à partir de ces figures que le contrôleur multi-boucle IMC proposé 
(IMC FOPI-FOF) atteint une bonne performance en terme du pourcentage de 
dépassement (%) et du temps de stabilisation (sec) par rapport aux performances de 
contrôleur multi-boucles (F-SDSP et SDSP) présenté par Chuong et al. [120, 75] et le 
contrôleur décentralisée fournie par Chananchai et al. [121]. 

En plus d'évaluer l'indice de performance en boucle fermée, l'intégrale de la valeur 
absolue de l'erreur (IAE) a été pris en compte, et il est observé que les performances 
avec le contrôleur proposé donnent de bons résultats et supérieurs par rapport aux les 
autres contrôleurs. 

Aussi à partir de la Figure 4-12 (a et b), il est clair que le contrôleur IMC fractionnaire 
proposé génère et donne une procédure de contrôle très fluide parmi ces méthodes et 
ceci est une confirmation des résultats obtenus. 

D'après l'analyse des performances du système dans le domaine fréquentielle les 
valeurs de la marge de phase et de fréquences de coupure obtenues sont très proches 


de celles spécifiées par le modèle de référence de chaque boucle. 


e En tant que contribution de notre part dans cette thèse et à partir des résultats obtenus 
par notre méthode proposée et comparaison des performances avec les derniers 
travaux atteints par les chercheurs, nous avons confirmé l'efficacité et le succès de la 
méthode proposée malgré la difficulté de mise en œuvre sur des systèmes 
multivariable. 

e _. Nous pouvons également dire que l'utilisation du contrôleur multi-boucle IMC 
proposé avec un découpleur inversé pour les systèmes MIMO (2x2) en utilisant la 
technique proposée offre de bonnes performances et robustesse. Cela est dû à 
l'utilisation de la fonction de transfert idéale de Bode en boucle fermée comme modèle 


de référence. 


4.5 Conclusion 


La nouveauté dans ce chapitre, est l'approche proposée pour résoudre le problème de 
commande fractionnaire des systèmes multivariables, par le choix d'une structure de 
commande par modèle interne (IMC) de ce type de systèmes, nous avons transformés un 
système MIMO (2x2) en deux systèmes SISO indépendants (deux boucles indépendants) en 
introduisant un découplage inversé pour résoudre le problème des interactions, puis à chaque 
système SISO indépendant a été réduit au modèle NIOPTD-I équivalent dans le but de la 
conception et la commande des système. La fonction de transfert de filtre f(s) égal à la 
fonction idéale de Bode en boucle fermée a été choisie comme modèle de référence lors de la 
synthèse du contrôleur. Les paramètres du contrôleur proposé ( IMC-FOPI-FOF controller) 
sont déterminés pour satisfaire aux spécifications de fréquence de coupure et marge de phase 
désiré pour chaque boucle. Les résultats de la simulation montrent clairement que les 
spécifications du domaine temporel et l'indice de performance du réponse du système de 
contrôle avec le contrôleur multi-boucles IMC d'ordre fractionnaire proposé donnent une 
bonne performance et acceptable en comparaison avec les autre contrôleurs (F-SDSP , C- 


SDSP et C-PID). 


Chapitre 5: Synthèse des contrôleurs PID fractionnaires pour 
les systèmes SISO/MIMO basée sur schéma du prédicteur de 
smith (SP). 


5.1 Introduction: 


Le temps de retard est l'un des obstacles qui conduisent à la déstabilisation du système de contrôle 
en boucle fermée ou, le contrôleur PID classique reste incapable d’imposer les performances désirés. 
Pour réduire les problèmes causés par le temps de retard dans le processus à commander, diverses 
stratégies de contrôle appropriées sont utilisées. Le meilleur d'entre eux est la stratégie de conception 
analytique de Prédicteur de Smith (SP) proposée en 1957 [122-123]. Il est également considéré aussi 
le schéma de contrôle le plus réussi et le plus largement utilisé pour les systèmes qui possèdent un 
temps de retard important. 


Le rôle principal de Prédicteur de Smith (SP) est qu'il renforce le système de contrôle et que le 
temps de retard est éliminé de l'équation caractéristique des systèmes en boucle fermée [122-123-124- 
125-126-1271]. 

Dans notre thèse, L’objectif de l’utilisation un nouveau modèle ordre fractionnaire de type II 


(NIOPTD-II) pour assurer une meilleure approximation des systèmes d'ordre élevé. Dans la première 
partie de ce chapitre, nous présentons une nouvelle structure de prédicteur de Smith (SP) proposée 
pour les systèmes multivariables avec découplage simplifié visant à éliminer le problème temps de 
retard qui caractérise ces systèmes. D'autre part, nous utilisons la fonction de transfert idéale de Bode 
en boucle fermée comme un modèle de référence et sa structure est équivalente à celle de SP. Dans la 
deuxième partie de ce chapitre, le contrôleur résultant est basé sur le principe de la synthèse directe. 


Appelé un: SP-FOF-FOPID contrôleur. 


5.2 Synthèse du contrôleur basé sur l'approche de prédicteur de Smith 
(SP) pour les systèmes (SISO) 


5.2.1 Principe de Prédicteur de Smith 


Dans la littérature, la méthode de prédicteur de smith vise à compensé l'effet des retards 
dans la boucle de régulation. où le temps de retard est théoriquement éliminé de la fonction de 
transfert en boucle fermée. Cette stratégie utilise un modèle de processus pour prédire le 
comportement de système et améliorer ainsi les performances du contrôleur. La structure 


fonctionnelle du prédicteur de Smith est donnée dans la figure 5.1. 


Figure 5-1:Structure fonctionnelle du prédicteur de Smith. 


où: y(s) : est la sortie de processus; u: est le signal de commande; R(S): signal de référence 
:d(s): est la perturbation d'entrée. 
, G(s) est le processus réel, G-(s) est le contrôleur d'ordre fractionnaire, G,(s)* est le modèle 
parfait, il peut être défini comme suit : 

Gu(s) = Gy(s)e”°"s (5.1) 
avec, Gm(s) est la partie libre de temps de retard du modèle, et em représente le temps de retard du 
modèle. 


La fonction de transfert en boucle fermée est donnée par 


FC) | Gc(s)GP(s) (5.2) 

R(S)  1+Gc(s)(Gmts) + Gp(s) — Gy(s)) 

Si le modèle est parfait c'est-à-dire Gp(s) = Gy(s)e-ÊmS Cela élimine la fonction du retard du 

dénominateur dans l'équation et se réduit a 
Y(S) _ Gc(S)Gm(s)e"?ms 
R(S)  1+Gc(S)Gu(S) 


Le schéma équivalent du prédicteur de Smith peut être représenté par figure 5.2. 


(5.3) 


Cp) 


Figure 5-2:Schéma équivalent du prédicteur de Smith. 


Le contrôleur peut être déterminé en utilisant l'approche de synthèse directe (DS) [17], ce qui donne : 


GS) = — (rs), 5.4 
Fe 7 Gm(s) e-ms — (RS) Le 
d 


Nous utilisons dans cette étude la fonction de transfert en boucle fermée proposée par Bode. 


Y(S)\ _ LB(S) os 1 _— 
Cor Dre Ro (6.5) 


où: Lg (s }est la fonction de transfert idéale de Bode en boucle ouverte et exprimée dans l'équation. 


(5.5). 
Le contrôleur G-(s) est obtenu en substituant (5.5) en (5.4). 
1 
REZ) “e 


La structure du système en boucle fermée désirée est illustrée à la figure 5.3. 


Figure 5-3:Schéma bloc de commande du système en boucle fermée désiré. 


5.2.2 Synthèse des contrôleurs SP-PID-FOF/ SP-FOPID-FOF 


5.2.2.1 Cas d'un système décrit par un modèle rationnel du premier ou deuxième ordre à 
retard (FOPDT/SOPDT) 
e Modèle du premier ordre plus un retard (FOPDT): 
On considère le modèle rationnel (d'ordre entier) du premier ordre plus un retard 


Gu(s)' donné par la fonction du transfert suivante : 


K 
Gus) = et (5.7) 
m 


En substituant l'équation (5.7) (sans délai) dans (5.6), on obtient le contrôleur Gç(s) sous la 


forme: 


TP GMCS)  Kmte TP Km 


Gc(s) = 


1 Ts +1 1 T 1 
" ( ) (5.8) 


TnS 
On note qu'il peut être mise sous forme d'un contrôleur PI d'ordre entier (PI) en série avec un filtre 


d'ordre fractionnaire (FOF), (appelé : contrôleur SP-PI-FOF) : 


1 
Ge(s) = F(s). (x (1 F —) (5.9) 
N 1 Tm 
ou: FC) = Er ke = T1 = Tom 


e Modèle du deuxième ordre plus un retard (FOPDT): 
On considère le modèle rationnel (d'ordre entier) du deuxième ordre plus retard, Gy(s)* 


donné par la fonction du transfert suivante : 


Km 
Ce = — 
nSY = psr1 


En substituant l'équation (5.10) (sans délai) à (5.6), on obtient le contrôleur Gç(s) à la forme: 


— Om S (5.10) 


7 1 _As?+Bs+1 1 (14245) ai 
Fe TB Gus)  KnteS  T. SB-1'\K Bs B° en 
C m 


On note qu'il peut être mise sous forme d'un contrôleur PID d'ordre entier (PID) en série avec un filtre 


d'ordre fractionnaire (FOF), (appelé : contrôleur SP-PID-FOF) : 


1 
Ge(s) = F(s). (x. (1 + —+ 5) (5.12) 
TS 
à 1 B A 
Où: FC) = ke tu = B,tn=>. 


5.2.2.2 Cas d'un système décrit par un modèle d’ordre non entier à retard de type I ou II 
(NIOPTD-I/NIOPTD-IT) 


e Modèle d’ordre non entier à retard de type I (NIOPTD-I): 
On considère le modèle d’ordre non entier à retard de type I, G,(s)* donné par la fonction du 


transfert suivante : 


* — Km — Om S 
Gu(s) TT,s+1 Lu 0<u<2; (5.13) 


En substituant l'équation (5.13) (sans délai) dans (5.6), on obtient le contrôleur Gç(s) à la 


forme: 


1 Tsk +1 1 T, 1 
_ (1+; ) 


G — = = 
c(s) TeSBGuts)  KmTeSP T SEK» 


. (5.14) 
m 


On note qu'il peut être mise sous forme d'un contrôleur PI d'ordre fractionnaire (FOPT) en série avec 


un filtre d'ordre fractionnaire (FOF), (appelé : contrôleur SP-FOPI-FOF) : 


1 
Gc(s) —= F(s). (x (1 + 5) (5. 15) 
à 1 Tr a 
ou: FC) = pr» Ke = T1 = Tom 


e Modèle d’ordre non entier à retard de type II (NIOPTD- ID): 
On considère le modèle d’ordre non entier à retard de type II, G/(s)* donné par la fonction 


du transfert suivante : 


* Km — Om S 
Gm(s) Lane 1<H2 <3; H>H: (5.16) 
En substituant l'équation (5.16) (sans délai) dans (5.6), on obtient le contrôleur Gç(s) sous la 


forme: 


Gc(s) = 


TP GUN) Kms Te PAT UKm 


1 As 2 + Bsu +1 1 B 
( Bs"1 


À u2- 
LE ci «)) (5.17) 


On note qu'il peut être mise sous forme d'un contrôleur PID d'ordre fractionnaire (FOPID) en série 


avec un filtre d'ordre fractionnaire (FOF), (appelé : contrôleur SP-FOPID-FOF) : 


Ces) = F(s). (x (1 + _ + m5") (5.18) 


N 1 B A 
ou: FC) = pa ke = Bt. 
Systèmes d'ordre supérieur 
(modèle original ) 
Première démarche Deuxième démarche Troisième démarche Quatrième démarche 


: + 


Approximation du modèle original par un | {Approximation du modèle original par un Approximation du modèle original par un | | Approximation du modèle original par un 
modèle d'ordre entièr à retard FOPDT [modéle d'ordre entierà retard SOPDT modèle fractionnaire de type NIOPTD-I| | modèle fractionnaire de typeNIOPTD-II 
\ 5 ni Y 

Calcul du contrôleur fractionnaire 
sur la base du modèle réduit 


Calcul du contrôleur fractionnaire 
sur la base du modèle réduit 
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Figure 5-4: Démarches à suivre pour la synthèse de contrôleurs d'ordre fractionnaires basé sur 


l'approche de prédicteur de Smith (SP) pour les systèmes d'ordre supérieur. 


5.3 Synthèse du contrôleur basé sur l'approche de prédicteur de Smith 
(SP) pour les systèmes MIMO: 


Dans cette partie, la structure de commande basée sur le prédicteur de Smith (SP) pour 
les systèmes SISO sera généralisée à la commande des systèmes multivariables (systèmes 
TITO) avec temps de retard afin d'améliorer les performances du système en boucle fermée. 
L'effet de toutes les boucles est considéré comme une perturbation sur la boucle de commande 
considérée, d'où l'idée d'utiliser découplage, pour séparer les interactions entre les variables 


du système MIMO. 


La méthode proposée est réalisée en combinaison avec des techniques de découplage qui 
simplifient la matrice de fonction de transfert des systèmes MIMO et abordent ensuite un 
schéma SP pour chaque boucle unique comme celui des systèmes SISO [17]. Le grand 
avantage de cette combinaison est qu'elle supprime complètement les termes de retard des 
éléments diagonaux de la matrice apparente découplée qui ont été approchées par des modèles 


d'ordre non entier plus un temps de retard de II type (NIOPTD-IT). La méthode proposée est 


basée sur l'approche de synthèse directe dans laquelle les contrôleurs PI/PID d'ordre 
fractionnaires sont conçus sur la base de la fonction de transfert idéale de Bode en boucle 


fermée désirées [16]. 


5.3.1 Prédicteur de Smith (SP) avec découplage simplifié. 


Dans ce travail, la structure de prédicteur de Smith avec découplage simplifié est utilisée 
pour traiter les principaux problèmes des processus multivariables, par exemple les 
interactions entre les variables de processus et les temps de retard multiples. Après la 
conception du découpleur simplifie, le procédé TITO est transformé en deux systèmes SISO 
indépendants. Étant donné que les boucles sont indépendantes dans ce cas, la structure du 
prédicteur de smith (SP) pour les systèmes SISO est appliquée avec les mêmes étapes pour 


chaque boucle indépendante. L'ensemble de la structure générale SP est illustré à la figure. 5.5. 


Figure 5-5:Structure générale du prédicteur de Smith pour le système TITO avec découplage 


simplifié. 


Le découpleur convertit les systèmes TITO en deux boucles individuelles, comme illustré à la 


Figure. 5.6. 
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Figure 5-6:Les boucles indépendantes résultantes 


Dans la Figure. 5.6, GA (s)",, et Gn(s)",, sont les modèles d'ordre fractionnaires, avec: 


Gn(SY'i 0 | 


: _[qu(s) 0 | 
Q(s) = G(s).D(s) = | … al — | 0 Gis) 


0 (5.19) 


où: Gnls)" is = q1ui(s) et GS)" 32 = q22(s). 
Les modèles d'ordre fractionnaires (les modèles de type II : NIOPD-IT) sont équivalents. 
Pour les éléments diagonaux ont été trouvés en suivant la méthodologie proposée à la section 


suivante. 


5.3.2 Modélisation des systèmes découplée 


Les systèmes découplée résultants sont des systèmes d'ordre élevé dans ce cas , il est difficile de 
concevoir directement les contrôleurs pour les éléments diagonaux de matrice Q(s) . Et donc, 
plusieurs techniques différentes ont été proposées pour approximation des éléments diagonaux 
(q11(S),q22(s)) par des modèles du premier ou deuxième ordre avec retard (modèles d'ordre entier : 
FOPDT/SOPDT) [7,95]. Mais parfois ces modèles ne sont pas exacts, fiables ou satisfaisants, d’où la 
nécessité d’un outil puissant pour contourner Les difficultés de modéliser les systèmes d’ordre élevé. 

On propose d’exploiter l’approche d'algorithme d'optimisation afin d’arriver à un modèle final 
capable de représenter convenablement les systèmes linéaires d’ordre élevé dans ses différents 
domaines de control. Par conséquent, les systèmes indépendantes sont approximées en modèles d'ordre 
fractionnaire équivalents (modèle d’ordre non entier à retard de type II (NIOPTD-ID [117, 16, 


133]). Ce modèle a été exprimée dans l'équation (5.20) ci-dessus 


K 
GO a su Bat I 


mii 


e-0miis où:1<H>;; <3; Hoi >Uii;i-12 (5.20) 


Les paramètres de chaque modèle sont donnés par le vecteur: 


X1 = [Km11 A11 B11 Om11 M211 Hi 11] pour le premier modèle Gm(s)';; 


(5.21) 
X2 = [Km22 A22 B22 0m22 M1 22 U2 22 | pour le deuxième modèle Gy(s)",, 
Dans l'étape de modélisation, les paramètres du modèle fractionnaires sont déterminés à partir 
de résoudre le problème d'optimisation proposés à l'aide de l'algorithme d'optimisation PSO. 
La fonction de cout est formulée en utilisant le critère du minimum de l'erreur 
quadratique moyenne (MSE) E,,;;, qui est donné par la valeur de l'écart entre les deux 
modèles du système réel Q(s) :;, et du modèle fractionnaire proposé Gn(S)ii. 


Jin (Km 11, À 11, B 11: Om 112 Hi 11 > 2 11 ) 


n 
1 
= => On — ÿ,)? pour le premier modèle . (5.22) 
RE 
Jmin (Ki 221 À 22, B 22, Om 22; M1 22 , U2 22 ) 
1 n 
= => On — ÿ,)? pour le deuxième modèle. (5.23) 
ni 
avec: 
Y1 = dus); Y2 = qa2(s); VE TIORE Ÿ2 = Gu(s) 


En outre, les paramètres du modèle fractionnaire proposés sont déterminés en résolvant le 
problème d'optimisation dans les conditions suivant : 


Kn min < Km 11,22 < Ka max 
0 < A4 1122 < Amax 
0 <B1122 < Bmax 
0 < On 11,22 < On max 
0 <Hij1122 <2 


1<Hz1122 <3;aVeC'Uzi > Mi 


(5.24) 


La structure de l'algorithme d'approximation est illustrée à la figure. 5.7 (pour plus de détails 


sur la modélisation du systèmes, voir [16]). 
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Figure 5-7:Le schéma fonctionnelle de la modélisation les systèmes découplée à modèles 
fractionnaires équivalents (NIOPTD-IT). 
5.3.3 Synthèse du contrôleur d'ordre fractionnaire (SP-FOPID-FOF): 


En utilisant la structure de contrôleur mentionnée ci-dessus (Figure. 5.7), les systèmes 
multivariables (systèmes de deux entrées et de deux sorties, en anglais: TITO system) 
deviennent des systèmes multi-boucles. Pour chaque boucle, un contrôleur correspondant doit 
être conçu pour répondre aux exigences de ses réponses en boucle fermée. Dans cette étude, 
une nouvelle structure de régulateur PID fractionnaire [128, 129]. est proposée pour chaque 
boucle, appelée SP-FOPID-FOF qui a été mentionnée ci-dessus l'équation (5.20). Ainsi, le 


contrôleur principal de chaque boucle est l'équation suivante : 


Gaits) = 


1 A;s ni2 + B;suit +1 1 B; 1 À pui 
a <Bi-uil' (1 + as + — sht HE ) (5. 25) 
ToiSPiGmi(s) Ki TaSP! Te SAM | Kou Bis? B; 


,avec:(i = 1; 2) 
La méthode de synthèse du contrôleur SP-FOPID-FOF se réalise en trois étapes : 
1-Analyse des interactions et choix de la configuration de commande : 
La matrice des gains relatifs (RGA) permet d'évaluer le niveau des interactions entre les 
variables du système multivariable et de choisir les couples entrée/sortie. 
2- Modélisation des systèmes découplée en modèles d'ordre fractionnaire équivalents de type 
IL (NIOPTD-ID). 
3-Chaque paramètre de contrôleur Prédicteur de Smith SP-FOPID-FOF de l'équation (5.25) 
est alors calculé comme suit : 
Les deux paramètres 7,; et Bi sont les paramètres de modèle de référence imposé à la boucle 


(ui, Yi). Cependant, dans notre cas, on propose d'imposer un comportement fractionnaire pour 


chaque boucle en choisissant la fonction idéale de Bode en boucle fermée comme modèle de 


référence [118, 119]. 


Y(s); Lg (Ss ); 1 
Cu) LRO = 1,2 (5:26) 
R(S); d 1 + Le (s ); 1 + TeiSPi 


Les deux paramètres 7,; et Bi sont choisis alors dans le domaine temporel en spécifiant la 
valeurs numériques du dépassement M, (%) ou bien le temps de réponse désirés t, (2%) pour 
chaque boucle, ceci en utilisant les relations de l'équations (1.38 et 1.39) dans la section (1.8) 


du première chapitre. 


5.4 Simulation et Application 


5.4.1 Application aux systèmes monovariable (SISO) 


Dans cette section, deux systèmes SISO réels plus un retard [130] seront choisis afin de valider 
les performances du contrôleur fractionnaire proposé (SP-FOPID-FOF) en fonction du modèle 
d'ordre fractionnaire de type II (NIOPDT-I) : 


e Premier cas : 


Cette fois ci, on appliquera la méthode de commande proposé de SP sur un système réel à 


décalage multiple SISO, défini dans [131], par la fonction de transfert suivante : 


GO = 
pour exécuter l'algorithme proposé dans la phase de modélisation, les paramètres initiaux sont 
sélectionnés et organisés dans le tableau suivant : 

Tableau 5-1:Les paramètres initiaux pour exécuter l'algorithme PSO dans la phase de 


modélisation (Cas d'un système monovariable): 


Paramètres Modèle fractionnaire de type II (NIOPDT-ID [132], 
Taille de la population 20 
nombre maximum d'itérations 120 
Nombre de particules 6 
Limite inférieure et limite —5<K, <5 
supérieure des variables 0 < A,B < 20 
0<06,, <20 


1<u2 <3;:0<h1 <2;U2 >: 


Coefficient cognitif Ci 1.5 
coefficient social C> 2 
Coefficient d'inertie W 1 


Le vecteur des valeurs finales des paramètres du modèle fractionnaire obtenu est : 


X= [Km À B 6» UH1 U2l 
=[1 10.674 15.6948 4.7771 2.1018 1.0296|, pour le modèle Gy(s)* 
Alors, le modèle fractionnaire de type II est le suivant: 
1 
10.674152 1018 + 5.69485102% + 1° 


Où l'écart de l'erreur obtenus dans la phase de modélisation entre système réelle et le modèle est 


Jmin = 1.3 x10 6. 


_ -47771s 


Gin (S)propose 


Notez que le modèle d'ordre entier ou dite un modèle de 2ème ordre avec retard (noté 
SOPTD 1[133]), proposé par Wang [131-134], pour le même système précédent, est d'défini 


par 
1 


EnSwnos = DE —5.39s 
9 (2.43s + 0.995)(2.43s + 0.995) 


e 


Les réponses indicielles du processus réel G,(s) et le modèle fractionnaire proposé sont 


illustrés à la figure. 5.8- 


Réponse indicielle 
T 


1.2 T T T T T I T T 
1 hé 
0.8 + 
2 
= L_ 
S 0.6 
an 
0.4 
Original 
—-—- Reduced-Order Model (NIOPTD-II)_ proposed by Boudjehem 
0.2 —-—- Reduced-Order Model (SOPTD)_proposed by Wang 
O L L l L L L L 1 L 
0 10 20 30 A0 50 60 70 80 90 100 
Temps (sec) 
Figure 5-8:Réponse indicielle : ---Processus réel, ---modèle d'ordre non entier parfait 


(NIOPTD-IT) proposé, --- modèle d'ordre entier parfait (SOPDT1[131 ]). 


Les valeurs désiré de M, (%) et t, (2%) sont 2.38% et 9.5s. respectivement. En utilisant 
équations (15-17). Nous obtenons les paramètres du filtre fractionnaire (7, et B) sont : T. = 
1.7609 ; P = 1.088. 

Maintenant, le contrôleur de prédicteur de Smith d'ordre fractionnaire (ou prédicteur de 
Smith FOF-FOPID ) peut être obtenu en utilisant la méthode de synthèse directe (DS) decrite 


dans la section 3, sous la forme: 


+1. s7aasto7)| 


1 
: Li 
- 6948 L 5.694851-0296 


Ge(s) = Loosovss { 


Les résultats obtenus sont comparés aux résultats donnés en [130]. 
Le contrôleur PID fractionnaire a été donné par Boudjehem et al [130]. Et il a la forme 


suivante: 


1.996 
= 0.5210 
C(s) = 0.2510 (1 +-0ose7 + 19.735 ) 


De l'autre coté, le contrôleur classique proposé par Wang dans [131] est 


C(s) = 0.46 (1 + _— +1. 235) 
et 
6.255 + 2.56 
"ss + 2.29) 


Le résultat de la réponse indicielle du système en boucle fermée commandé par le contrôleur SP-FOF- 


C(s) = 


FOPID et ceux des techniques de Boudjehem et de Wang sont illustrés à la figure. 5.9- 
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Figure 5-9:Réponses indicielles du Premier système contrôlé (système SISO) avec le 


contrôleur proposé et celles données dans [130]. 


Les meilleures valeurs des indices de performance obtenues par le contrôleur SP-d'ordre fractionnaire 
(prédicteur de Smith FOF-FOPID) sont comparés avec ceux donnés par Boudjehem et de Wang [130- 
131] dans le Tableau 5.2. 


120 


Tableau 5-2:L'analyse les performances de la réponse indicielle par le contrôleur proposé et 


ceux donnés dans [130-131]: 


Méthode de réglage Temps de réponse | Dépassement (%) IAE ISE 
(secondes) M; 
Ts (2%) 
Contrôleur SP-FOF- 14.69 2.98 3.141 0.9244 
FOPID proposé 
D. boudjehem et al. 15.11 2.25 4.575 1.485 
[130] Contrôleur SP- 
FOPID 
Wang et al. [131] 23.32 2.98 3.609 1.373 
mismatched SP-PID 
Wang et al. [131] 18.88 3.068 17.93 10.07 
perfectly matched SP 


Selon la figure.5.9. et le tableau 5.2., le contrôleur fractionnaire proposée donne des performances 
mieux que celles obtenus par le contrôleur FOPID proposé par boudjehem [130] et le contrôleur PID 
classique proposé par Wang [131] en termes de tous les indices de performance (IAE ,ISE et ITAE). 
D'après les données du Tableau (5.2), les valeurs du dépassement M, (%) et du temps de réponse 
t,(2%) obtenues sont très proches de celles spécifiées par le modèle de référence. 
Nous remarquons qu'il y a un meilleur rejet des perturbations de la réponse indicielle avec le 
contrôleur proposé contrairement aux autres contrôleurs comme le montre la figure 5.9. 

e Deuxième cas: 


Dans ce deuxième cas, la fonction de transfert avec retard pur du système réel défini dans [131-134], 


est: 
1 
G = —4s 
5-55 
Dans la phase de modélisation, le modèle d'ordre fractionnaire de type II obtenu est: 
1.0454 


_ —5.273 


où l'écart de l'erreur obtenus dans la phase de modélisation entre système réelle et le modèle est 
Jmin = 8.1539 x 1074 
Le modèle d'ordre entier SOPDT proposé par Wang dans [131-134] est défini par : 


1 
G —_ _ __— 
m(S)wang © 645 + 0.999 ° 


Figure.5.10.llustre les réponses indicielles du processus réel G,(s) et le modèle d'ordre 


—5.79s 
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fractionnaire. 
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Figure 5-10:Les réponses indicielles : ---Processus réel, ---modèle d'ordre réduit (NIOPTD- 


ID) proposé, --- modèle d'ordre réduit (SOPDT!{131]) proposé par Wang. 


Dans cet exemple, les valeurs souhaitées pour M, (%) et t, (2%) sont 2.0248% et 5.988s. 
successivement. En utilisant les équations (15-17), on obtient les paramètres du filtre 
fractionnaire (t,. et B) comme suit :T. =1.1609 ;  B = 1.078. 

La fonction de transfert du contrôleur d'ordre fractionnaire SP-FOF-FOPID obtenu est sous 


la forme : 


Gc(s) = (3 6681 (1 + +1. sagstesse)) 


1 
1. 1609 50-078 ° 3.668150-?7008 
Les résultats obtenus sont comparés aux résultats [130]. 


D'autre part, le contrôleur PID classique a été donné par Wang et al [130] , est sous la forme 


suivante 
C(s) = 1.12 (1 + 3775 + 0.94) 
et 
5.03s + 3.06 
CCS) = ——— 
s(s + 2.48) 


La figure 5.11 montre la réponse indicielle du système commandé par le contrôleur d'ordre 


fractionnaire proposé et ceux des techniques de Wang et al [130]. 
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Figure 5-11:Réponses indicielles du deuxième système contrôlé (système SISO) avec le 


contrôleur proposé et celles données dans [130]. 


On remarque d'après la figure 5.11, que le contrôleur fractionnaire proposé à la capacité d'atteindre 
l'amélioration souhaitée par rapport aux autres contrôleurs Ceci est traduit par l'obtention d'une 
réponse temporelle caractérisée par un dépassement faible MP(%) et un temps de réponse t,(2%) 
réduit. Ils sont aussi très proches de celles spécifiées par le modèle de référence. 

À partir des résultats de simulation obtenus grâce à notre stratégie proposé de SP comme contribution 
de notre part à ce travail en le comparant avec les derniers travaux atteints par les chercheurs, nous 
avons confirmé l'efficacité et le succès de la méthode proposée aux systèmes SISO. Nous pouvons 


maintenant l'étendre pour inclure des systèmes MIMO 2x2 qui ont des difficultés à mettre en œuvre. 


5.4.2 Application aux systèmes multivariable (MIMO) 


Le modèle mathématique de la colonne de distillation est un système multivariable de deux entrées 
et de deux sorties 2x2 (A modified Alatiqi subsystem distillation column [6]), est exprimé par la 


matrice de transfert suivante: 


—0.51 1.68 
—_———— ,#L"75s LS 
to (32s + 1)2(2s + 1) (28s + 1)2(2s + 1) 
LE 1.25 2.8s 4.78 1155 


A8s+1(6s+10° 
La matrice des gains relatifs est obtenue à l'aide de l'équation (2.10) comme suit: 


7.2167 —6.2167 
—6.2167 7.2167 


En utilisant l'équation (2.20), on obtient l'indice de Niederlinski : NI = 0,1386. 


(43.65 + DS +0 
RGA(G(0)) = | 


D'après l'analyse des éléments de la RGA et comme NI est positif, la configuration de commande 


appropriée est : la sortie y. est contrôlée par le signal de commande u; et la sortie y, est contrôlée par 
le signal de commande u. 

D'après la section 2 « De la partie : découplage des systèmes multivariable», la matrice de transfert du 
découpleur simplifié peut être donnée par: 


3.294(32s + 1)? 


1 
(28s + 1)? 


D(S)2x2 = 


0.262 (48s + 1)(55 + 1)e 


es 
(43.65 + 1)(9s + 1) 


Les éléments diagonaux résultants (les systèmes découplées) de H(s) = G(s)2x2 X D(s)2x2 = 


diag(h::1{s), h22(s)) sont : 


Ge —0.51e 75 : 0.439(48s + 1)(55 + 1)e 3655 
US G25+1)C5+1)  (285+ 1)2 (43.65 + 1)(05+ 1)(25 + 1) 
4.78e 6655 4.118(32s + 1)2e 25 
h22(s) = 


(485+1)(55+1) (28s+1)2(43.6s + 1)(95+1) 
Dans la phase de modélisation des systèmes découplées, les paramètres initiaux sont 


sélectionnés et organisés dans le tableau suivant : 


Tableau 5-3:Les paramètres initiaux pour exécuter l'algorithme proposé dans la phase de 


modélisation (Cas d'un système multivariable): 


Paramètres Le modèle fractionnaire G,, 11 CS) pour Le modèle fractionnaire 
le premiere sous-système découplé Gm 22 (S) pour le deuxième sous- 
système découplé 
Taille de la population 20 
nombre maximum d'itérations 100 
Nombre de particules 12 
Limite inférieure et limite supérieure des —5 <Ky11 <5 —5<K/y22 <5 
variables 0 < A;1,B11 < 100 0 < 4,2, B22 < 100 
0 < 0,511 < 100 0 < 0,22 < 100 
1<u: <3; 0<u:; <2; 1<u: <3; 0<u; <2; 
H2 > HU: H2 > HU: 
Coefficient cognitif Ci 1.5 
coefficient social C2 2 
Coefficient d'inertie W 1 


Les modèles d'ordre fractionnaire de type IT équivalents aux systèmes découplées ont été 


obtenus en suivant la méthodologie proposée dans la section 3 et donnée comme suit: 


—0.071 


G (s) 1, -803059s 
mA proposed 34, 08775169 + 80.4603510754 + 1 


e 


100 


0.6624 
En 22(Sproposed © EG 054S T7 + 76.06855108 + 1° 


Dans cette étape de modélisation, l'écart de l'erreur obtenu entre les deux systèmes réels 


—634377s 


découplés et les modèles d'ordre non entier NIOPTD-II est la suivante : 


Jmin = 9.29 x107% pour le premier modèle et /min = 1.119 x 10! pour le deuxième 
modèle. 

Les modèles du premier ordre avec retard (noté: FOPTD) équivalents aux systèmes découplés 
obtenus par la méthode de tavakoli [95], sont donnés par : 


—0.071 


Gm 11 (S)tavakoli = 58.761521 —77.24s 
0.662 
Gm 22 (S)tavakoli = 77.24s11 —61981s 


Les réponses indicielles du modèle d'ordre réduit de type NIOPTD-II proposé sont données sur la 


Figure 5.12, pour évaluer l'affinité avec réponse de chaque sous-système découplé. 
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Figure 5-12:La réponse indicielle des systèmes découplés original et leurs modèles d'ordre 


fractionnaire équivalents NIOPTD-II. 
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On peut voir que les réponses indicielles du modèle NIOPTD-II sont très proches de leurs systèmes 
découplées équivalents. Ceci est considéré comme une indication de l'efficacité de la stratégie de 
l'approximation proposée. 

Par ailleurs, le contrôleur PI décentralisé d'ordre entier proposé par Tavakoli et al. ([95]) basé sur 
l'utilisation de la méthode de réglage NDT pour les modèles FOPDT réduit, est donné par : 


0.092 
—6.393 — Fu 0 
C(S)vpr = .012 
0 0.499 + 


De l'autre coté, le contrôleur PID décentralisé à d'ordre entier proposé par Wang’s dans [6] est donné 
par : 


9 
—2.176 — + 0.807s 0 


C(S)wang S 


5 
0 0.154 + — 0.1175s 


Les valeurs de M,,(%) et t; (2%) imposes sont 6,1372% et 450s. En utilisant la méthode de synthèse 


directe (DS) décrite dans la section 3, la fonction de transfert du contrôleur de prédicteur de Smith est 


donné . 
C(s) 
1 
sise (-1183.2487 (+ areas + 2428750000) . 
f nn (usa 837 (1 a 0.65805°71)) 
121.291950122 76.0685s1053 


Les réponses indicielles du système en boucle fermée contrôlé par le contrôleur SP-fractionnaire et 
les contrôleurs PI/PID classiques sont illustrés dans la Figure. 5.13 (bouclel (Figure. 5.13.a) et 
boucle? (Figure. 5.13.b)). 
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Réponse indicielle 
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Figure 5-13:Réponses indicielles du système contrôlé (système 2x2 ) avec le contrôleur 


proposé et ceux donnés dans [95]. 


Les valeurs des indices de performance ISE et IAE et Les spécifications dans le domaine temporel ( 


(M, et T; (2%)) pour les sorties contrôlées basées sur trois techniques de réglage différentes sont 


reportées dans le tableau 5.4. 


Tableau 5-4:Indices de performance et spécifications de domaine temporel: 


Méthode de réglage Entrée (u)- | Temps de réponse Dépassement 2 DE 
Sortie (y) (secondes) T, (2%) M, (%) 
Contrôleur SP-FOF- TEA 445.20 13 101.2 | 48.11 
FOPID proposé [17] U_Y2 404.88 6.52 95.23 54.31 
Tavakoli et al. [95] Ui_Y1 SEL 1110 250.7 175.9 
NDT-PI W_y2 63306 LUXE 2627 169 
Wang et al. [6] Ui_Y1 7465 2.73 3585 227.7 
Contrôleur PID U_Y2 685.8 1.09 4075 2443 


e À partir de la Figure. 5.13 et le tableau 5.4., nous constatons que le contrôleur SP-FOF- 


FOPID conçu par la méthode proposée permet de donner des meilleures performances en 


termes de pourcentage de dépassement et de temps de réponse par rapport aux contrôleurs 
PI/PID classiques. 


e On observe également que les valeurs des indices de performance (IAE et ISE) avec le 
contrôleur d'ordre fractionnaire donnent des résultats très supérieurs par rapport au contrôleur 
classique (C-PT). 

e Il ressort de cette figure que le contrôleur proposé est plus puissant, et qu'il n'est pas affecté 
par les perturbations. 

e On peut dire que l'utilisation d'un contrôleur SP-FOF-FOPID proposé basé sur le schéma de 
prédicteur de Smith pour les systèmes MIMO avec un découpleur simplifié utilisant une 


approche de synthèse directe offre de bonnes performances et robustesse. 


Pour démontrer l'efficacité de la stratégie proposé et atteindre les spécifications dans le domaine 
temporel les plus exigeantes, nous avons spécifié plusieurs valeurs différentes désiré du 
M, (%) et T; (2%). 

Proposed( Smith Predictor FOF-FOPID })} 


—— —  y1 (PO=S.51 : ts (2%)-450) 
——— y1 (PO=9.79 ; ts (2%)-600) 
———y1 (PO=15.45 ; ts (2%)}-700) 
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(a) 
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———» y2 (PO=9.79 : ts (2%)=600) 
——. y2 (PO=15.45 ; ts (2%)=700) 
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Figure 5-14:Réponses indicielles en boucle fermée du système contrôlé avec différentes 


valeurs désiré de M, (%) et T; (2%). 


104 
es 


Le tableau suivant montre les paramètres des contrôleurs obtenus selon les valeurs 


imposées au modèle de référence. 


Tableau 5-5:Les paramètres des contrôleurs proposés en fonction des valeurs imposées du 


modèle de référence. 


T, (2%) = 4505 


1 
H 80.460351-0754 


+ 0.4237s061%) 


Les valeurs imposées - , 1 
Contrôleur SP-FOF-FOPID proposé de type: G-(s) = F(s).| k. (1 +—+ Tps) 
pour M,(%) et Us 
T, (2%) Boucle 1 Boucle 2 
ie ou : 1133.2437 (1 : 114 837(1 
117.501750-0867 * ‘ 117.501750-1091 ° . 


1 
di 76.068551-053 


+ 0. 65950715) 


M, = 9.79% 


T, (2%) = 6005 


1 
230.52489 50.166 * 


( 1133.2437 (1 


1 
L 80.460351-0754 


+ 0.4237s061%) 


230.5248501890 * 


1 
(araes7(1 


1 
bi 76.068551-053 


+ 0. 65900715) 


M, = 15.45% 


T,; (2%) = 700s 


1 
121.2919502526 


( 1133.2437 (1 


1 
. 80.4603510754 


+ 0.4237s061%) 


121.29195 02750” 


1 
(aiaes7(1 


1 
EL 76.0685s1-053 


+ 0. 65908075) 


“D'après les réponses indicielles (Figure. 5.14.a et b), nous voyons que si les valeurs désiré de 
M,(%) et T; (2%) sont petites, il en résulte de petites valeurs de tç et B , ce qui conduit à la 
réponse plus rapide du système et donne un bon rejet de la perturbations. 

“ Nous mentionnons également que les valeurs de filtre d'ordre fractionnaire (ou de la première 

partie du contrôleur proposé) que nous obtenons correspondent aux valeurs de M,(%) et 

T; (2%) imposées comme suit : 

e Si M, = 9,79% et T; (2%) = 600s imposé de cette manière, les valeurs de filtre d'ordre 
fractionnaire obtenues sont T7. = 231,0137 et ff = 1.242. 

e Si M, = 15,45% et T; (2%) = 700s imposés comme tels, on obtient les valeurs du 
fractionnaires de filtre Tr. —468.8246 et f —1,328, tout en gardant la même seconde partie 
(contrôleur PID fractionnaire) du régulateur proposé inchangée dans tous les cas car ses 
paramètres sont déterminés à partir du modèle d'ordre fractionnaire (NIOPDT-IT). 

e Il est clairement démontré que cette efficacité de commande et ces améliorations ont été 


obtenues en satisfaisant les performances désirées aux systèmes multivariables. 
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e L'efficacité de l'utilisation du contrôleur proposé est évidente dans le rejet de la 
perturbation et de l'atténuation des interactions entre les entrées/sorties du système 


contrôlé. 


5.5 Conclusion 


Dans ce chapitre, nous avons étudié deux problèmes d’une classe des systèmes multivariables avec 
retard. La méthode de commande basé sur le prédicteur de smith avec découplage simplifié 
développée dans ce chapitre est capable de renforcer le système de contrôle et l'élimination des effets 
de temps de retard et d'interférences entre les entrées et les sorties du système. 

Les résultats de simulation ont montré que la stratégie de synthèse du contrôleur proposé appliquée 
au système monovariable SISO ainsi qu'au système complexe multivariable MIMO a prouvé son 
efficacité et sa robustesse par rapport aux différentes techniques de conception. Alors que les 
contrôleurs d'ordre fractionnaire basés sur la structure de prédicteur de Smith appliqués aux modèles 
d'ordre non entier ont contribué à obtenir les résultats les plus efficaces en terme de pourcentage de 


dépassement M,(%) et de temps de Réponse ( t, ) pour ceux introduits par d’autre méthodes. 


NN 
[e] 
O 


Conclusion générale: 


L'objectif principal des travaux de recherche que nous avons menés dans le cadre de cette 
thèse est la contribution à la synthèse des lois de commande à retour fractionnaire pour les 
systèmes multivariables carrés. 

La première contribution de ce travail est de considérer le couplage fort qui existe entre 
les variables d'entrée et de sortie des systèmes multivariables comme des perturbations dans 
les boucles de régulation, et donc. Nous avons réduit l'influence des interactions en 
introduisant des techniques de découplages et en obtenant des boucles indépendantes. Dans la 
seconde contribution, nous nous sommes appuyés sur le développement de certaines 
méthodes de synthèse pour les contrôleurs d'ordre fractionnaires issus des systèmes SISO et 
leur adaptation aux systèmes TITO , en plus de proposer de nouvelles techniques de synthèse 
plus robuste et plus efficace, notamment dans le cas des systèmes MIMO stables. Les résultats 
obtenus de nos recherches sont organisés selon l'ordre des chapitres. 

Nous avons présentés les notions de base des calcul fractionnaire, en particulier les 
opérateurs d’ordre fractionnaire et leurs techniques d'approximation, et Nous avons fait deux 
exemples de simulation de opérateurs fractionnaires approchés par la représentation diffusive. 
Aussi, nous avons introduit quelques notions sur le principe de la méthode CRONE basée sur 
la fonction idéale de Bode, et étudié le comportement de la fonction en boucle fermée 
résultante dans les domaines temporel et fréquentiel afin de montrer ses caractéristiques et 
l'avantage de l'utiliser comme un modèle de référence dans la synthèse du contrôleur d'ordre 
fractionnaire. 

Nous avons présentées différent méthodes des synthèses des contrôleurs classiques 
trouvées dans la littérature, avec des simulations faites pour voir les limites de efficacité ces 
techniques. 

Nous avons appliqué la synthèse optimal basé sur l'algorithme d'optimisation au le 
modèle réel de colonne de distillation est un système de deux entrées et deux sorties (2x2) 
avec découplage simplifié. Nous avons également appliqué seconde simulation à un système 
multivariable à trois entrées et trois sorties (3X3). Enfin, nous avons conclu que l'algorithme 
PSO est très efficace pour optimiser les paramètres des contrôleurs PI d'ordre fractionnaire 
basé sur l'approximation diffusive, ce qui assure un excellent comportement du système 
multivariables en comparaison des contrôleurs PI classique (NDT et SIMC). Les résultats de 
simulation ont montré que la synthèse optimal peut être implémenté d’une façon plus efficace aux les 


systèmes MIMO par l’approche proposée. 
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Nous avons développé une stratégie de commande pour les modèles d'ordre non entier basé 
sur l'équivalence entre la structure de commande par modèle interne IMC et la structure classique à 
retour unitaire, nous avons choisi la fonction de transfert de filtre de causalité est égal à la fonction 
idéale de Bode en boucle fermée comme un modèle de référence. Par la suite, Le nouveau contrôleur 
résultant peut être décomposé en un contrôleur PID d'ordre fractionnaire cascadé avec un filtre 
fractionnaire. Dans la partie simulation, nous fournissons un exemple illustratif d'un système 
monovariable pour la stratégie de synthèse proposée, Puis nous avons étendu cette technique pour 
l'adapter à la première fois aux les systèmes multivariables et cela a été fait en deux étapes: le 
premièr, nous avons séparés le système MIMO en utilisant une technique de découplage inversée qui 
produit deux systèmes indépendants d'ordre élevé, puis sont approximé avec des modèles d'ordre 
fractionnaires de type I, le deuxième étape, nous avons synthétisé le contrôleur IMC-FOF-FOPI multi- 
boucle d'ordre fractionnaire pour les boucle indépendants. Les résultats de la simulation ont montré 
que la technique de contrôle proposée est plus efficace et acceptable par rapport aux autres techniques 
de synthèse trouvées dans la littérature. 
Afin d'améliorer les performances du schéma de prédicteur de Smith (SP) en boucle fermée pour les 
systèmes à temps mort, nous avons proposé model d'ordre non entier de type II (NIOPTD-ID) . Et eu 
l'idée d'utiliser la fonction de transfert idéale de Bode en boucle fermée comme un modèle de 
référence et sa structure est équivalente à celle de SP. Par la suite, la structure de commande basée sur 
le prédicteur de Smith (SP) pour les systèmes SISO a été développé pour contrôler les systèmes 
multivariables avec des retards, Dans ce contexte, nous avons proposé une nouvelle structure de 
controle bien adaptée à cette but. Les résultats de la simulation ont montré que la stratégie de synthèse 
direct du contrôleur proposé a prouvé son efficacité et sa robustesse par rapport aux différentes 


méthodes de réglage existantes malgré la difficulté de contrôle du système en boucle fermée. 
Les perspectives pouvant être envisagés, parmi les quels on peut citer : 


- Extension de l’implémentation des stratégies de commande décentralisée d'ordre fractionnaire 
développées dans cette thèse aux les systèmes multivariables (de plus de deux entrée et deux sorties ) 
et les systèmes multivariables non carrée . 

- J'étendrai également cette méthodologie aux systèmes MIMO fractionnaires et chercherai des 
moyens appropriés pour découplés ceux systèmes. 

- Le développement d’autres méthodes de synthèse des contrôleurs fractionnaire et plus adaptées aux 
systèmes MIMO instables. 

- Contribution à la synthèse des lois de commande centralisés d'ordre fractionnaires pour les 


systèmes multivariables. 
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